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Abstract 

We establish the long time soliton asymptotics for the translation invariant nonlinear 
system consisting of the Klein-Gordon equation coupled to a charged relativistic particle. 
The coupled system has a six dimensional invariant manifold of the soliton solutions. We 
show that in the large time approximation any finite energy solution, with the initial state 
close to the solitary manifold, is a sum of a soliton and a dispersive wave which is a solution 
of the free Klein-Gordon equation. It is assumed that the charge density satisfies the Wiener 
condition which is a version of the "Fermi Golden Rule" . The proof is based on an extension of 
the general strategy introduced by Soffer and Weinstein, Buslaev and Perelman, and others: 
symplectic projection in Hilbert space onto the solitary manifold, modulation equations for 
the parameters of the projection, and decay of the transversal component. 
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1 Introduction 



Our paper concerns the problem of nonlinear field-particle interaction. A charged particle ra- 
diates a field which acts back on the particle. This interaction is responsible for some crucial 
features of the dynamics: asymptotically uniform motion and stability against small perturba- 
tions of the particle, increase of the particle's mass and others (see [TJ [TTJ EHj)- The problem 
has many different appearances: a classical particle coupled to a scalar or Maxwell field, and 
coupled Maxwell-Schrodinger or Maxwell-Dirac equations, general translation invariant nonlinear 
hyperbolic PDEs. In all cases the goal is to reveal the distinguished role of the soliton solutions, 
i.e. traveling wave solutions of finite energy. Let us note that the existence of the soliton solutions 
is proved for nonlinear Klein-Gordon equations with a general nonlinear term j^Ej, and for the 
coupled Maxwell-Dirac equations [T3]. 

One of the main goals of a mathematical investigation is to study soliton type asymptotics 
and asymptotic stability of soliton solutions to the equations. First results in this direction have 
been discovered for the KdV equation and other completely integrable equations. For the KdV 
equation, any solution with sufficiently smooth and rapidly decaying initial data converges to a 
finite sum of soliton solutions moving to the right, and a dispersive wave moving to the left. A 
complete survey and proofs can be found in [T2] . 

For nonintegrable equations, the long time convergence of the solution to a soliton part 
and dispersive wave was obtained first by Soffer and Weinstein in the context of ?7(l)-invariant 
Schrodinger equation [3IJ . The extension to translation invariant equations was obtained 

by Buslaev and Perelman [SJ E] for the ID Schrodinger equation, and by Miller, Pego and Wein- 
stein for the ID modified KdV and RLW equations, j23 EE I2H] • The techniques introduced by 
Weinstein jH] play a fundamental role in the proofs of all these results. 

In 01] the long time convergence is obtained for the ID translation invariant and £7(1)- 
invariant nonlinear Schrodinger equation. It is shown there that the following asymptotics hold 
for any finite-energy solution ip(x,t) with initial data close to a soliton ip Vo (x — v^t — ao)e lu3ot : 



Here the first term on the right hand side is a soliton with parameters v±, a±, u>± close to t>o, ao, 
ujq, the function Wo(t)ip± is a dispersive wave which is a solution to the free Schrodinger equation, 
and the remainder r± (x, t) converges to zero in the global L 2 -norm. Recently Cuccagna extended 
the asymptotics (jl.l|) to nD Schrodinger equations with n > 3, PI ITU]. 

We establish the asymptotics similar to (jl.lj) for a scalar real- valued Klein-Gordon field ip(x) 
in H 3 coupled to a relativistic particle with position q and momentum p governed by 

tp(x, t) = 7r(x, t), 7r(x, t) = A^(x, t) — m 2 ip(x, t) — p(x — q(t)), x E IR 3 , 



where m > (the case m = is degenerate and will be considered elsewhere). This is a Hamilto- 
nian system with the Hamiltonian functional 



ip(x,t) = ip v± (x - v±t - a±)e luJ±t + W (t)ip± + r±(x,t), t -> ±oo. 



(1.1) 




(1.2) 
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The first two equations for the fields are equivalent to the Klein-Gordon equation with the source 
p(x — q). The form of the last two equations in (jl.2|) is determined by the choice of the relativistic 
kinetic energy ^/l + p 2 in (jl.3|) . Nevertheless, the system (jl.2|) is not relativistic invariant. 

We have set the maximal speed of the particle equal to one, which is the speed of wave 
propagation. This is in agreement with the principles of special relativity. Let us also note that the 
first two equations of (jl.2|) admit the soliton solutions of finite energy, ip v (x — vt — a), 7i v (x — vt — a), 
if and only if \v\ < 1. 

The case of a point particle corresponds to p(x) = 5(x) and then the interaction term in the 
Hamiltonian is simply ip{q). However, in this case the Hamiltonian is unbounded from below 
which leads to the ill-posedness of the problem, also known as ultraviolet divergence. Therefore 
we smooth the coupling by the function p(x) following the "extended electron" strategy proposed 
by M. Abraham pQ for charges coupled to the Maxwell field. In analogy to the Maxwell-Lorentz 
equations we call p the "charge distribution" . Let us write the system (jl.2j) as 

Y(t) = F(Y(t)), t G 1R, (1.4) 

where Y{t) := (tp(x,t),n(x,t),q(t),p(t)) (below we always deal with column vectors but often 
write them as row vectors). The system (jl.2|) is translation- invariant and admits the soliton 
solutions 

Y a ^ v (t) = (ip v (x — vt — a),7t v (x — vt — a),vt + a,p v ), p v = v / V 1 — v 2 (1.5) 

for all a, v G IR 3 with \v\ < 1 (see (|2.7J) . ()2.10|) ). where the functions ip v , ir v decay exponentially 
for m > (main difficulty of the case m = is provided by very slow decay of the functions). The 
states S a>v := 1^,(0) form the solitary manifold 

S := {S a>v :a,ve IR 3 , \v\ < 1}. (1.6) 

Our main result is the soliton-type asymptotics of type (jl.lj) for t — > ±oo, 

(ifj(x,t),7r(x,t)) ~ (ip v± (x - v±t - a±),ir v± (x - v±t - a±)) + W (*)*± (1-7) 

for solutions to (jl.2j) with initial data close to the solitary manifold S. Here Wo(t) is the dynamical 
group of the free Klein-Gordon equation, are the corresponding asymptotic scattering states, 
and the remainder converges to zero in the global energy norm, i.e. in the norm of the Sobolev 
space if 1 (]R 3 ) © L 2 (IR 3 ). For the particle trajectory we prove that 

q(t)^v±, q(t)~ v± t + a±. (1.8) 

The results are established under the following conditions on the charge distribution: p is a real 
valued function of the Sobolev class if 2 (IR 3 ), compactly supported, and spherically symmetric, 
i.e. 

p, Vp, VVp G L 2 (JR 3 ), p(x) = for |x| > R p , p{x) = Pl {\x\). (1.9) 

We require that all "modes" of the wave field are coupled to the particle, which is formalized by 
the Wiener condition 

p(k) = (2vr)- 3 / 2 / e lkx p{x)dx ^ for all k G IR 3 . (1.10) 
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It is an analogue of the "Fermi Golden Rule" [Tj |HJ E01 1201 EH 100] : the coupling term p(x — q) 
is not orthogonal to the eigenf unctions e tkx of the continuous spectrum of the linear part of the 
equation. As we will see, the Wiener condition (jl.lUj) is very essential for our asymptotic analysis 
(see Remark 115. 5J1 . Generic examples of the coupling function p satisfying 1)1.9)1 and (Jl.lUj) are 
given in [21]. 

Remark 1.1 Physically, the Wiener condition means the strong coupling of the particle to the 
field which leads to radiation of the particle. This radiation results in the relaxation of the 
acceleration q(t) — > 0, t — > ±oo which provides the asymptotics 1)1.7)1 and 1)1.8)1 . Note that the 
soliton solutions do not radiate, and the radiation of the particle manifests itself in the decay of 
the deviation of the solution from the solitary manifold (see ()1.18|) below). 

The problem under investigation was studied earlier in the following two different situations 
A and B: 

A. The asymptotics 

Q(t)->v±, (i/j(x,t),7r(x,t)) ~ (il) v ±(x-q(t)),n v ±(x-q(t))) (l.H) 

were proved in [20] in the case m = 0, under the Wiener condition (|1.10j) . for all finite energy 
solutions, without the assumption that the initial data are close to S. This means that the 
solitary manifold is a global attractor for the equations (|1.2|) . However, the asymptotics (|l.llj) 
were established only in local energy semi-norms centered at the particle position q(t). This 
means that the remainder in 1)1.11)1 may contain a dispersive term, similar to the middle term in 
the right hand side of (jl.7|) . whose energy radiates to infinity as t — > ±oo but does not converge 
to zero. A similar result is established in JO] fo r coupled Maxwell-Lorentz equations. 

B. The asymptotics 1)1.110 . and an analogue of the asymptotics 1)1.7)1 in the global energy norm, 
were established in [TH] (resp., [TO]) also for all finite energy solutions, in the case m = (resp., 
m > 0), under the smallness condition on the coupling function, HjoHxamss <C 1. The similar 
results are established in [T71 [01] (resp., [19, ) for the coupled Maxwell-Lorentz equations with 
a moving (resp., rotating) charge. Let us stress that the asymptotics p. 8)1 for the position was 
missing in the previous work. 

Let us comment on main difficulties in proving the asymptotic stability of the invariant man- 
ifold S and justifying ([1.70 . 1)1.8)1 . The method of [101 [20] is based on the Wiener Tauberian 
Theorem, hence cannot provide a rate of convergence in the velocity asymptotics of (Jl. 8)1 which is 
needed to prove (jl. 7)1 and the position asymptotics of 1)1.8)1 . Also the methods of [TB*) IT71 ITHl IT9"] 
are applicable only for a small coupling function p(x), and do not provide the position asymptotics 
in (OJl . 

Our approach develops a general strategy introduced in [01 El I2H1 120] for proving the asymptotic 
stability of the invariant solitary manifold S. The strategy originates from the techniques in [UJ 
and their developments in [Oil [021100] m the context of the [/(l)-invariant Schrodinger equation. 
The approach uses the symplectic geometry methods for the Hamiltonian systems in Hilbert spaces 
and spectral theory of nonself adjoint operators. 

The invariant manifolds arise automatically for equations with a symmetry Lie group |H EH Kl- 
in particular, our system 1)1.2)1 is invariant under translations in 1R 3 . The asymptotic stability of 
the solitary manifold is studied by a linearization of the dynamics (jl. 4)1 . The linearization will 
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be made along a special curve on the solitary manifold, S(t), which is the symplectic orthogonal 
projection of the solution. Then the linearized equation reads 

X(t) = A(t)X(t), t G IR, (1.12) 

where the operator A(t) corresponds to the linearization at the soliton S(t). Furthermore, we 
consider the "frozen" linearized equation ()1.12j) with A(ti) instead of A(t). The operator A(t\) 
has zero eigenvalue, and the frozen linearized equation admits secular solutions linear in t (see 
(|6.24Jl ). The existence of these runaway solutions prohibits the direct application of the Liapunov 
strategy and is responsible for the instability of the nonlinear dynamics along the manifold S. One 
crucial observation is that the linearized equation is stable in the symplectic orthogonal complement 
to the tangent space T$. The complement is invariant under the linearized dynamics since the 
linearized dynamics is Hamiltonian and leaves the symplectic structure invariant. 



Our proofs are based on a suitable extension of the methods in jSJ 12E1 EH] ■ Let us comment 
on the main steps. 

I. First, we construct the symplectic orthogonal projection S(t) = UY(t) of the trajectory 
Y(t) onto the solitary manifold S. This means that S(t) G S, and the complement vector Z(t) : = 
Y{t) — S{t) is symplectic orthogonal to the tangent space Ts(t) for every t G IR: 

z(t)\r S{th tern. (i.i3) 

So, we get the splitting Y{t) = S(t) + Z(t) and we linearize the dynamics in the transversal 
component Z(t) along the trajectory. 

The soliton component S(t) = Sb(t), v (t) satisfies a modulation equation. Namely, in the parametriza- 

tion £(i) = (c(t),v(t)) with c(t) := b{t) — I v(s)ds, we have 

Jo 

i{t)=N 1 (t(t),Z{t)), \Ni(t(t),Z{t))\<C\\Z(t)\\ 2 _ fi , (1.14) 

where || • stands for an appropriate weighted Sobolev norm. 

On the other hand, the transversal component satisfies the transversal equation 

Z(t) = A(t)Z(t) + N 2 (S(t), Z(t)), (1.15) 

where A(t) = Asa), and iV^S^t), Z(t)) is a nonlinear part: 

||iV 2 (^),^))||^<C||Z(i)|| 2 _ /3 , (1.16) 

where || ■ \\p is defined similarly to || ■ Let us note that the bound (|1.16|) is not a direct 

consequence of the linearization, since the function S(t) generally is not a solution of ()1.4|) . The 
modulation equation and the bound ()1.14|) play a crucial role in the proof of f)l. 16j) . 

II. The linearized dynamics (|1.12|) is nonautonomous. First, let us fix t = t\ in A(t) and 
consider the corresponding "frozen" linear autonomous equation with A(t\) instead of A(t). We 
prove the decay 

\\X(t)\\-fi ^ ^+|^ » teM ( L1? ) 
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of the solutions X(t) to the frozen equation for any X(0) G Zs 1 where S\ := S(tx), and Zs 1 is the 
space of vectors X which are symplectic orthogonal to the tangent space T Sl . Let us stress that 
the decay holds only for the solutions symplectic orthogonal to the tangent space. Basically, the 
reason of the decay is the fact that the spectrum of the generator A(tx) restricted to the space 
Zg 1 is purely continuous. 

III. We combine the decay f|1.17|) with the bound ()1.14|) through the nonlinear equation ((1.15)1 . 
This gives the time decay of the transversal component 

W>ll-,< f™ ten, 

if the norm ||Z(0)||/3 is sufficiently small. One of the main difficulties in proving the decay (jl.l8|) is 
the non-autonomous character of the linear part of ((1.15)1 . We deduce the decay from the equation 
p. 15)1 written in the "frozen" form 

Z{t)=A{t 1 )Z{t) + [A{t)-A{t 1 ))Z{t) + N 2 {S{t),Z{t)), <t<h, (1.19) 

with arbitrary large t± > 0. 

IV. The decay ()1. 18)1 implies the soliton asymptotics ((1.7)1 and ((1.8)1 by the known techniques 
of scattering theory. 

Remarks 1.2 i) The asymptotic stability of the solitary manifold S is caused by the radiation of 
energy to infinity which appears as the local energy decay for the transversal component, ((1.18)1 . 

ii) The asymptotics () 1.7)1 can be interpreted as the collision of the incident soliton, with a trajectory 
w_t+a_, with an incident wave Wo(t)^-, which results in an outgoing soliton with a new trajectory 
v + t + a+, and a new outgoing wave Wo{t)^f + . The collision process can be represented by the 
diagram of Fig. 1. It suggests to introduce the (nonlinear) scattering operator 

S: ( V _,a_,*_)^(t; + ,a + ,* + ). (1.20) 

However, the domain of the operator is an open problem as well as the question on its asymptotic 
completeness (i.e. on its range). 



Remarks 1.3 i) The strategy of (51 El I2E I2H] was further developed in the papers [3 El 13 
Uni 123 OH HH3 EH] ■ Let us stress that these papers contain several assumptions on the discrete 
and continuous spectrum of the linearized problem. In our case a complete investigation of the 
spectrum of the linearized problem is given under the Wiener condition and there is no need for 
any a priori spectral assumptions. 

ii) Note that the Wiener condition is indispensable for our proof of the decay p. 17)1 . but only in 
the proof of Lemma 115.31 Otherwise we use only the fact that the coupling function p(x) is not 
identically zero. The other assumptions on p can be weakened: the spherical symmetry is not 
necessary, and one can assume also that p belongs to a weighted Sobolev space rather than having 
a compact support. 

Our paper is organized as follows. In Section 2, we formulate the main result. In Section 3, we 
introduce the symplectic projection onto the solitary manifold. The linearized equation is defined 
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Figure 1: Wave - particle scattering. 



in section 4. In Section 6, we split the dynamics in two components: along the solitary manifold 
and in the transversal directions, and we justify the estimate (|1.14j) concerning the tangential 
component. The time decay of the transversal component is established in sections 7-10 under 
an assumption on the time decay of the linearized dynamics. In Section 11, we prove the main 
result. Sections 12-18 fill the gap concerning the time decay of the linearized dynamics. In 
Appendices A and B we collect some routine calculations. 

Acknowledgments The authors thank V.Buslaev for numerous lectures on his results and fruitful 
discussions, and E.Kopylova for reading the paper and making many useful remarks. 

2 Main Results 

2.1 Existence of Dynamics 

To formulate our results precisely, we need some definitions. We introduce a suitable phase space 
for the Cauchy problem corresponding to (jl.2j) and (jl.HJ) . Let H° = L 2 be the real Hilbert space 
L 2 (IR 3 ) with scalar product (-, ■) and norm || ■ \\ L 2, and let H 1 be the Sobolev space H 1 = {tp 6 
L 2 : |V^| £ i 2 } with the norm HV'IIh 1 = HV^Hl 2 + HV'IU 2 - Let us introduce also the weighted 
Sobolev spaces H^, s — 0, 1, a G IR with the norms H^lUa : = ||(1 + M ) "Cliff's • 

Definition 2.1 i) The phase space S is the real Hilbert space H 1 © L 2 © H 3 © IR 3 of states 
Y = [ijj, Ti, q,p) with the finite norm 

\\y\\ £ = IMIffi + 1Mb + M + bl- 

ii) S a is the space Hi © Hi © IR 3 © IR 3 with the norm 

\\Y\\a = II *1k = IM|l.a + |M|0,a + \q\ + \ P \. (2.1) 
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Hi) T is the space H 1 © L 2 of fields F = (ip,7r) with the finite norm 

\\F\\r = IMU 1 + IklU 2 - 
Similarly, T a is the space H\ © H® with the norm 

\\F\\a= H^lk = IM|l,a+Ho,a- (2-2) 

Note that we use the same notation for the norms in the space T a as in the space S a defined 
in (|2.ip . We hope it will not create misunderstandings since T a is equivalent to the subspace 
of £ a which consists of elements of S a with zero vector components : q = p = 0. It will be 
always clear from the context if we deal with fields only, and therefore with the space or with 
fields-particles, and therefore with elements of the space £ a . 

We consider the Cauchy problem for the Hamilton system (jl.2}) which we write as 

Y(t) = F(Y(t)), teJR: Y(0) = Y . (2.3) 

Here Y(t) = (tl>(t),ir(t),q(t),p(t)), Y = (ipo,7i ,q ,p ), and all derivatives are understood in the 
sense of distributions. 

Proposition 2.2 [15J Let / TO) hold. Then 

(i) For every Y e £, the Cauchy problem \2. S\) has a unique solution Y(t) G C(H, £). 

(ii) For every t G E, the map U(t) : Yq t— > Y(t) is continuous on £. 
(Hi) The energy is conserved, i.e. 

H(Y(t))=H(Y ), ten, (2.4) 

and the velocity is bounded, 

|g(t)|<«<l, teM, (2.5) 

where v = v(Y ) . 

The proof is based on a priori estimates provided by the fact that the Hamilton functional (|1.3|) 
is bounded from below. The latter follows from the bounds 

h < ^UWh + <iM- - v)) < ^^mh + h\p\\b, (2-6) 



2m 2UrnL ~ 2 lir11 ^ ' Nr ' rv ^" ~ 2 " rili ' 2 
which imply also that 8 is the space of finite energy states. 

2.2 Solitary Manifold and Main Result 

Let us compute the solitons (jl.5|) . The substitution to (jl.2|) gives the following stationary equa- 
tions, 

-v-Vip v (y) = w v (y) , -v-V-K v {y) = Aif) v (y) - m 2 ijj v (y) - p(y) 



Pv o = -J vMy)p(y)dy 



VT+p. 



(2.7) 
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Then the first two equations imply 

M>M ■= [-A + m 2 + (v V) 2 ]^(y) = -p(y), 



y e IR 3 . 



(2.8) 



For \v\ < 1 the operator A is an isomorphism if 4 (IR 3 ) — > H 2 (JR 3 ). Hence (jl.9|) implies that 

Mv) = -A~V(y) e ^ 4 (iR 3 )- (2-9) 

If v is given and \v\ < 1, then p v can be found from the third equation of (|2.7jl . Further, functions 
p and ■?/>„ are even by (|1.9jl . Thus, V?/v is odd and the last equation of ()2.7|) holds. Hence, the 
soliton solution (jl.5|) exists and defined uniquely for any couple (a,v) with \v\ < 1. 

The function can be computed by the Fourier transform. The soliton is given by the 
formulas 



i> v (x) 



7T v (x) 



J_ 

4:71 



i{y- x )w + (v-x)±\ 

V 



-v-Vij) v (x), Vv = lV 



(2.10) 



Here we set 7 = l/vl — v 2 and x — x» + x±, where xn \\v and x±-Lv for x E IR 3 . 
Let us denote by V := {v E H 3 : |u| < 1}. 

Definition 2.3 A soliton state is S(a) := (ip v (x — b),ir v (x — b),b,p v ), where a := (b,v) with 
b E IR 3 and dGV, 



Obviously, the soliton solution admits the representation S(o~(t)), where 

a(t) = (b(t),v(t)) = (vt + a,v). 



(2.11) 



Definition 2.4 A solitary manifold is the set S := {S(a) : a E E := IR 3 x V}. 

The main result of our paper is the following theorem. 

Theorem 2.5 Let hl.9p and 1 1.1 Up hold. Let /? > 3/2 and Y(t) be the solution to the Cauchy 
problem h2. with the initial state Yq which is sufficiently close to the solitary manifold: 

d :=dist £/3 (y ,«S) < I- (2.12) 
Then the asymptotics hold for t — > ±00, 

q(t) = v ± + 0(\t\- 2 ), q(t) = v±t + a ± + 0{\t\~^ (2-13) 

(ip(x,t),7r(x,t)) = {ip v± {x - v±t-a±),7T v± (x - v±t - a±)) + W (t)*± + r±(x,t) (2.14) 

with 

||r ± (t)||^ = 0(|t|- 1 / 2 ). (2.15) 
It suffices to prove the asymptotics (|2.14j) . 1)2.13)1 for t — > +00 since system (jl.2j) is time reversible. 



3 Symplectic Projection 

3.1 Symplectic Structure and Hamilton Form 

The system fjl.2|) reads as the Hamilton system 



Y = JVH(Y), J 



( 





1 





o\ 




-1 






















1 


\ 








-1 





Y =(V,7r, q,p)e£, (3.1) 



where T>7i is the Frechet derivative of the Hamilton functional ()1.3|) . Let us identify the tangent 
space of £, at every point, with the space £. Consider the symplectic form fl defined on £ by the 
rule 

£1 = dip(x) A d7r(x) dx + dq A dp. 



In other words, 

n(Y 1 ,Y 2 ) = (Y 1 ,JY 2 ), Y u Y 2 e£, (3.2) 

where 

(Yi, Y 2 ) := (ipuifa) + (vri, tt 2 ) + q x q 2 + pxp 2 
and (ipi,ip 2 ) = J ipi(x)ip 2 (x)dx etc. It is clear that the form £7 is non-degenerate, i.e. 

fi(Yi, Y 2 ) = for every Y 2 £ £ Y 1 = 0. 



Definition 3.1 T/ie symbol Yi \ Y 2 means that Y 1 E £,Y 2 E £, and Y 1 is symplectic orthogonal 
toY 2 , i.e. Q(Yi,Y" 2 ) = 0. 

ii) A projection operator P : £ — > £ is said to be symplectic orthogonal ifY\ \ Y 2 for Y\ G Kerf 
and Y 2 G ImP. 

3.2 Symplectic Projection onto Solitary Manifold 

Let us consider the tangent space Ts(a)S of the manifold S at a point S(o~). The vectors Tj : = 
d aj S(a), where d a . := and <9 CTj . + 3 := d v . with j = 1, 2, 3, form a basis in T a S. In detail, 



Tj = Tj (v) := d h .S(a) = ( -djip v (y) , -dj-K v (y) , ej , 
r i+3 = 7"j+3(^) := 9„ .^(o-) = ( d v ip v (y), d v ir v (y) , 0, d v p v 



j = 1,2,3, (3.3) 



where y := x — b is the "moving frame coordinate", e\ = (1,0,0) etc. Let us stress that the 
functions Tj are always regarded as functions of y rather than those of x. 
Formulas (|2.10|) and conditions (J1.9)) imply that 

Tj{v) G £ a , veV, j = 1, . . . , 6, Va G IR. (3.4) 

Lemma 3.2 The matrix with the elements Q(ri(v),Tj(v)) is non-degenerate for any v G V. 
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Proof The elements are computed in Appendix A. As the result, the matrix Q(ri,Tj) has the 
form 

n(«) := (0(71, Tj)),, , () = ( °_ n+(v) ^ +{V) ) , (3.5) 

where the 3 x 3-matrix Q + (v ) equals 

n+(v) =K + (1- v 2 )" 1/2 E + (1 - v 2 )~ 3/2 v <g> u. (3.6) 
Here X is a symmetric 3 x 3-matrix with the elements 

k 2 + m 2 + 3(kv) 2 f 2 A; 2 + m 2 + 3(kv) 



KlJ = J dk\Uk)\ hk 3 k2 + m2 _ {kv)2 = J dk\p(k)\ k ^ {k 2 + m 2_ {kv y r ( 3 - 7 ) 

where the "hat" stands for the Fourier transform (cf. (jl.lOjl ). The matrix K is the integral of the 
symmetric nonnegative definite matrix k®k = (kikj) with a positive weight. Hence, the matrix K 
is also nonnegative definite. Since the identity matrix E is positive definite and the matrix v ®v is 
nonnegative definite, the matrix Q + (v) is symmetric and positive definite, hence non-degenerate. 
Then the matrix Q(ti, Tj) is also non-degenerate. ■ 
Now we show that in a small neighborhood of the soliton manifold S a "symplectic orthogonal 
projection" onto S is well-defined. Let us introduce the translations T a : (ip(-),ir(-),q,p) \— > 
(ip(- — a),7r(- — a),q + a,p), a G IR 3 . Note that the manifold S is invariant with respect to the 



translations. Let us denote by v(p) := p/y 1 + p 2 for p e IR . 

Definition 3.3 i) For any a e IR and v < 1 denote by £ a (v) = {Y = (ip,ir,q,p) G £ a '■ \v{p)\ < 
v}. We set£(v) :=£ (v). 

ii) For any v < 1 denote by E(-u) = {a — (b, v) : b G IR 3 , |f | < ?)}. 

Lemma 3.4 Lei (Qj^ /ioW, a G IR andv < 1. T/ien 

i/iere exists a neighborhood O a (S) of S in £ a and a mapping II : O a («S) — > 5 such that II is 
uniformly continuous on O a (S) fl £ a (?J) in the metric of £ a , 

UY = Y for YeS, and Y-S\T S S, where S = UY. (3.8) 
ii) O a (S) is invariant with respect to the translations T a , and 

UT a Y = T a UY, for Y G O a {S) and a G IR 3 . (3.9) 
Hi) For any v < 1 there exists a v < 1 s.t. T1Y = S(a) with a G £({;) for Y G O a (S) fl £ a (v). 
iv) For any v < 1 there exists an r a {v) > s.t. S(a) + Z G O a (S) if a G £(£>) and ||^|| Q < r a (-y). 
Proof We have to find o = o~{Y) such that S(a) = UY and 

Q(Y - S(a), d a .S(a)) = 0, j = 1, . . . , 6. (3.10) 

Let us fix an arbitrary a G X and note that the system (J3.10|) involves only 6 smooth scalar 
functions of Y. Then for Y close to S(o~°), the existence of a follows by the standard finite 
dimensional implicit function theorem if we show that the 6x6 Jacobian matrix with elements 
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Mij (Y) = d ai Q(Y - S(a°),d aj S(a )) is non-degenerate at Y — S(a°). First note that all the 
derivatives exist by ([3.4)1 . The non-degeneracy holds by Lemma 13.21 and the definition (J3.3)) 
since Mij(S(o~ )) = — Q,(d ai S(a°), d a ,S(a )). Thus, there exists some neighborhood O a (S(a )) of 
S(a°) where II is well defined and satisfies (|3.8|) . and the same is true in the union 0' a (S) = 
U (J o eE C Q (5(a )). The identity (Q holds for y,T a y G 0«(<S), since the form and the manifold 
S are invariant with respect to the translations. 

It remains to modify 0' a (S) by the translations: we set O a (S) = U^^zT^O'^S) . Then the 
second statement obviously holds. 

The last two statements and the uniform continuity in the first statement follow by translation 
invariance and compactness arguments. ■ 

We refer to II as symplectic orthogonal projection onto S. 

Corollary 3.5 The condition implies that Y = S + Z where S = S(a ) = HY , and 

\\Z \\p<l. (3.11) 

Proof Lemma 13.41 implies that HY = S is well defined for small do > 0. Furthermore, the 
condition ()2.12|) means that there exists a point Si G S such that ||lo — Lillys = do- Hence, 
Y ,Si G Op(S) n S^(v) with some v < 1 which does not depend on d for sufficiently small d . 
On the other hand, HSi = Si, hence the uniform continuity of the mapping II implies that 
||<Si - S\\f3 — > as d -> 0. Therefore, finally, ||Zo||/3 = ||^o - S\\p < \\Y - SiWp + \\Si - S\\/3 < 
do + o(l) <^ 1 for small do. ■ 



4 Linearization on the Solitary Manifold 

Let us consider a solution to the system ([1.2)1 . and split it as the sum 

Y(t) = S(a(t)) + Z(t), (4.1) 

where a(t) = (b(t),v(t)) G S is an arbitrary smooth function of t G IR. In detail, denote Y = 
(if), Ti, q,p) and Z = ( 1 i r , II, Q, P). Then ()4.1j) means that 

if)(x,t) = if) v{t) (x-b(t)) + ^(x-b(t),t), q(t) = b(t) + Q(t) 

7r(x,t) = n v(t) (x-b(t))+Tl(x-b(t),t), p(t) = p v{t) +P(t) { ' } 

Let us substitute (J4.2)) to ([1.2)1 . and linearize the equations in Z. Below we shall choose S(a(t)) = 
IIF(t), i.e. Zit) is symplectic orthogonal to Ts(<j{t))S- However, this orthogonality condition is 
not needed for the formal process of linearization. The orthogonality condition will be important 
in Section 6, where we derive "modulation equations" for the parameters cr(i). 

Let us proceed to linearization. Setting y = x — b(t) which is the "moving frame coordinate", 
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(4.3) 



we obtain from (|4.2j) and ()1.2|) that 

i> = v-V v My)-b-VMy) + y(y,t)-b-\7V(y,t) = n v (y) + Tl(y,t) 

ii = v-V v ir v {y)-b-Vir v {y) + fl{y,t)-b-Vn{y,t) 

= Aip v (y) - m 2 ij v (y) + A* (j/, t) - m 2 ty(y, t) - p(y - Q) 

q = b + Q = - = 

y/1 + (p, + P) 2 

p = v-V vPv + P = -(V(il)v(y) + *(v,t)),p(y-Q)) 

The equations are linear in \l/ and II, hence it remains extract the terms linear in Q and P. First 
note that p(y - Q) = p(y) - Q ■ Vp(y) - N 2 (Q), where -N 2 (Q) = p(y - Q) - p(y) + Q ■ Vp{y). 
The condition (jl.9|) implies that for N 2 (Q) the bound holds, 



\N 2 (Q)\\o,f3 < Cfs(Q)Q 2 



(4.4) 



uniformly in \Q\ < Q for any fixed Q, where (3 is the parameter in Theorem 12.51 Second, the 
Taylor expansion gives 



Pv + P 



v/i + ( Pv + py 



v + v(P-v(v- P)) + N 3 (v,P), 



where v : = (1 + pi) 1 ^ 2 = a/1 — v 2 , and 



\N 3 (v,P)\<C(v)P 2 



(4.5) 



uniformly with respect to |i>| < v < 1. Using the equations (|2.7|) . we obtain from ()4.3|) the 
following equations for the components of the vector Z(t): 



(4.6) 



*(y, t) = U(y, t)+l- *) + (&-«)■ - « • V„Vi,(l/) 

n(y, t) = AV(y, t)-m 2 ^(y, t)+b • Vr%, t) +Q-Vp{y) + (b-v)-Vn v (y) -v ■V v n v (y) + N 2 
Q(t) =v{E-v® v)P + (v-b)+N 3 

P(t) = (*(y, *), Vp(y)) + (VMv), Q ■ Vp(y)) - « • V„p„ + N 4 (v, Z) 

where N A (v,Z) = (Vtp v , N 2 (Q)} + (W, Q • Vp) + (W,iV 2 (Q)}. Clearly, N±(y,Z) satisfies the 
following estimate 

|W 4 MI < Cp{p,v,Q)[Q 2 + 11*111,-/9^1], (4.7) 
uniformly in \v\ < v and \Q\ < Q for any fixed v < 1. We can write the equations ()4.6j) as 

= A(*)Z(t) + T(t) + JV(t), t e JR. (4.8) 
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Here the operator A(t) = A V)W depends on two parameters, v = v(t), and w = b(t) and can be 
written in the form 

/ ^ \ / w V 



A, 



n 

Q 



A — m 2 w ■ V 





Vp- 








T 



V 



/ 



N(a,Z) 



J 





(*\ 




n 




Q 






Z) are g 







N 2 (Z) 



(4.9) 



N 3 (v,Z) 
\ N A (v,Z) J 



(4.10) 



/ (w - v) ■ Vip v - v ■ V v ip v \ 
(w — v ) • X7tt v — V ■ V V 7T V 
v — w 
-v ■ V v p v 

where v = v(t), w = w(t), a = ait) = (b(t),v(t)), and Z = Z(t). Estimates (jO|) . (1431) and (gZZD 
imply that 

\\N{o,Z)\\ p <C{v,Q)\\Z\\%, (4.11) 
uniformly in a G S(-u) and H^H-/? < r -/?(^) f° r an y fixed t? < 1. 

Remarks 4.1 i) The term A(t)Z(t) in the right hand side of equation (j4.8|) is linear in Z(t), and 
N(t) is a /iz^/i order term in On the other hand, T(t) is a zero order term which does not 

vanish at Z(t) = since S(a(t)) generally is not a soliton solution if (|2.11j) fails to hold (though 
S(a(t)) belongs to the solitary manifold), 
ii) Formulas ()3.3j) and ()4.1U|) imply: 

3 

T W = - J2^ W - ^ + ^ Ti +3] ( 4 - 12 ) 
1=1 

and hence T(t) G Tgu^S, t G IR. This fact suggests an unstable character of the nonlinear 
dynamics along the solitary manifold. 



5 Linearized Equation 

Here we collect some Hamiltonian and spectral properties of the operator (|4.9jl . First, let us 
consider the linear equation 

X(t) = A VtW X(t), t G IR (5.1) 

with arbitrary fixed v G V = {v G IR 3 : \v\ < 1} and w G IR 3 . Let us define the space 
£+ := H 2 (WL 3 ) © H^WL 3 ) © IR 3 © IR 3 . 

Lemma 5.1 i) For any v G V and w G IR 3 equation \5.1\) can be represented as the Hamiltonian 
system (cf. iTH)) ). 

X(t) = JDH v , w {X(t)), t G IR, (5.2) 
where DTl VjV] is the Frechet derivative of the Hamiltonian functional 

1 



nr+ Wr + m i w| 



hw • Vtfdy + / p(y)Q • Wdy 



+ ip.^P-^g-V^^^-Vp^)), A- = (*,n,Q,P)6£. (5.3) 
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ii) The energy conservation law holds for the solutions X(t) G C 1 (IR, £ + ), 

H V;W (X(t)) = const, t G JR. 
Hi) The skew-symmetry relation holds, 

Q[A v>w Xi, X2) = —ft(Xi,A VtW X2), Xi,X 2 G £. 

Proof i) The equation (j5.1j) reads as follows, 

\ 



(5.4) 



(5.5) 





( * ^ 




d 


n 




di 


Q 






\p J 





( n + w-w 

A\l> - m 2 ^ + w-VU + Q-Vp 



B V P 



(5.6) 



The first three equations correspond to the Hamilton form since 

U + w -V^ = D n H v>w , A^-m^ + w VIl + Q- Vp = -DyH v , w , B V P = V P H ViW . 
Let us check that the last equation has also the Hamilton form, i.e. — (V\l/, p) + (Vip v , Q ■ Vp) 



—VqH V)W . First we note that — (dj^,p) 



pQ ■ V^/dx. It remains to show that 



(d^ v , Q ■ Vp) = d Q .-(Q ■ V^, Q ■ Vp). 



(5.7) 



Indeed, 



^(d j ip v ,Q-Vp) + -(Q-Vip v ,d j p) 
-{d^ v ,Q-Vp) + ^{d^ v ,Q-Vp) 



(5.8) 



where we have integrated twice by parts. Then ()5.7|) follows. 

ii) The energy conservation law follows by ()5.2)1 and the chain rule for the Frechet derivatives: 



jH v , w (X(t)) = (DH v , w (X(t)),X(t)) = (DH v , w (X(t)),JDH v , w (X(t))} = 0, 



t G 1R (5.9) 



since the operator J is skew-symmetric by (|3.1|) . and DTC V;W (X(t)) G £ for X{t) G £ + . 

iii) The skew-symmetry holds since A VjW X = JDH VjW (X), and the linear operator X \— > DH V)W (X) 

is symmetric as the Frechet derivative of a quadratic form. ■ 

Remark 5.2 One can obtain ()5.3|) by expanding 7i(Sb tV + X) to a power series in X up to second 
order terms. As a result, 7i VtW (X) is the quadratic part of the Taylor series complemented by the 
second integral on the right hand side of ()5.3|) arising from the left hand side of (|3.1|) . 

Lemma 5.3 The operator A v<w acts on the tangent vectors Tj(v) to the solitary manifold as fol- 
lows, 



A V)W [Tj(v)} = (w - v) -VTj(v), A VtW [r j+3 (v)] = (w - v) ■ Vr j+3 (v) +r j (v), j = 1,2,3. (5.10) 
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Proof In detail, we have to show that 



A, 



( -djip v \ 
-djir v 



\ / 



/ (v - w) ■ V djtp v \ 

(v — w) • Vdj7T v 





A v - 



( d Vj tp v \ 

9 Vj 7T v 



V d VjPv j 



( (w-v)- Vd v ijj v \ 



\ 



(w - v) ■ Vd Vj -K v 





/ 

/ -Ojlfjy \ 



J 



(5.11) 



V o ) 



Indeed, differentiate the equations (|2.7|) in bj and Vj, and obtain that the derivatives of soliton 
state in parameters satisfy the following equations, 



— v ■ Vdjip v = djir v , 
-djip v - v ■ Vd v .ij) v = d v .7r v 



—v ■ Vdj7r v = Adji() v — m 2 djip v — djp 



-djir v - v ■ Vd Vi -n v = Ad v ,ip v - m 2 d v .i/j v 



d Vi p v 



v 



3 



e ^ 1 -^- 1/2 + -( 1 _, 2) 3/ 2 



o = -(va.^,p) 



(5.12) 



for j = 1, 2, 3. Then (jo~TT| follows from (!5~T2l by definition of A in (jOJ). ■ 

We shall apply Lemma f5 . 1 1 mainly to the operator A v<a corresponding to w — v. In that case 
the linearized equation has the following additional specific features. 

Lemma 5.4 Let us assume that w = v G V . Then 

i) The tangent vectors Tj(v) with j = 1,2,3 are eigenvectors, and Tj +S (v) are root vectors of the 
operator A vv , that correspond to the zero eigenvalue, i.e. 



A v>v [Tj(v)) = 0, A V:V [r j+3 (v)] = Tj{v), j = 1,2,3. 
ii) The Hamiltonian function \5.§) is nonnegative definite since 



1 



\U + v W| 2 + |A 1/2 * - A~ 1/2 Q • Vpp 



dx + ^P- B V P > 0. 



(5.13) 



(5.14) 



Here A is the operator which is symmetric and nonnegative definite in L 2 (IR 3 ) for \v\ < 1, 
and A 1//2 is the nonnegative definite square root defined in the Fourier representation. 

Proof The first statement follows from f)5.10j) with w = v. In order to prove ii) we rewrite the 
integral in (j5.14j) as follows: 

-(II + v ■ W, n + v ■ W) + ^(A 1/2 ^ - A- 1/2 Q • Vp, A 1/2 vp - A~ 1/2 Q • Vp) 

_ _ 



i(A*, vi/) - (vi/, g • vp) + \(^Q ■ Vp, Q ■ Vp) 



(5.15) 
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since the operator A 1 / 2 is symmetric in L 2 (IR 3 ). Now all the terms of the expression (J5.15)) can 
be identified with the corresponding terms in ()5.3|) since 

^(Atf,tf) = ^([-A + m 2 + (v V) 2 ]^>, A-V = -V»> (5-16) 

by (EU) and (EH- ■ 

Remark 5.5 In Section 14 we will apply Lemma f5. 41 ii) together with energy conservation (|5.4|) 
to prove the analyticity of the resolvent (A v ^ v — A) -1 for Re A > 0. 

Remark 5.6 For a soliton solution of the system ()1.2|) we have b = v, v = 0, and hence T(t) = 0. 
Thus, equation 1)5.1)) is the linearization of system ()1.2j) on a soliton solution. In fact, we linearize 
(ll.2j) on a trajectory S(a(t)), where a(t) is nonlinear with respect to t, rather than on a soliton 
solution. We shall show below that T(t) is quadratic in Z{t) if we choose S(a(t)) to be the 
symplectic orthogonal projection of Y(t). In this case, (15. lj) is a linearization of (jl.2j) again. 



6 Symplectic Decomposition of the Dynamics 

Here we decompose the dynamics in two components: along the manifold S and in transversal 
directions. The equation (|4.8|) is obtained without any assumption on a{t) in 1)4.1)1 . We are going 
to specify S(a(t)) := YlY(t). However, in this case we must known that 

Y(t) G O a (S), t G H, (6.1) 

with some O a {S) defined in Lemma [3 .41 It is true for t = by our main assumption ()2.12)) with 
sufficiently small d > 0. Then S(a(0)) = UY(0) and Z(0) = Y(0) - S(a(0)) are well defined. 
We shall prove below that ()6.ip holds with a = —f3 if do is sufficiently small. First, the a priori 
estimate ()2.5)) together with Lemma [3.41 iii) imply that UY(t) = S(a(t)) with a(t) = (b(t),v(t)), 
and 

\v(t)\<v<l, teJR (6.2) 

if Y(t) G 0-p{S). Denote by r^p(v) the positive number in Lemma 13.41 iv) which corresponds 
to a = —f3. Then S(a) + Z G 0-p(S) if a = (b,v) with |u| < v and H^H-/? < r~p(v). Note 
that ()2.5j) implies ||Z(0)||_/3 < r-p(v) if d is sufficiently small. Therefore, S(a(t)) = HY(t) and 
Z(t) = Y(t) — S(a(t)) are well defined for small times t > 0, such that < r-p(v). This 

argument can be formalized by the following standard definition. 

Definition 6.1 Let t* be the "exit time", 

U = sup{t > : ||Z(s)||-/3 < r-p(v), 0<s<t}. (6.3) 

One of our main goals is to prove that £* = oo if do is sufficiently small. This would follow if 
we shall show that 

\\Z(t)\\-p<r- P (v)/2, 0<t<U. (6.4) 

Note that _ 

\Q{t)\ <Q:= r-p(v), 0<t<U. (6.5) 

Now by 1)4.11)1 . the term N(t) in ()4.8j) satisfies the following estimate, 

ll^(*)ll/J<C/9(«)ll^(*)ll-/3. 0<t<U. (6.6) 
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6.1 Longitudinal Dynamics: Modulation Equations 

From now on we fix the decomposition Y(t) = S(a(t)) + Z(t) for < t < t* by setting S(a(t)) = 
HY(t) which is equivalent to the symplectic orthogonality condition of type 1)3.8)1 . 

Z(t)\Tsw))S, 0<t<U. (6.7) 

This enables us drastically simplify the asymptotic analysis of the dynamical equation (J4.8)) for the 
transversal component Z(t). As the first step, we derive the longitudinal dynamics, i.e. find the 
"modulation equations" for the parameters cr(i). Let us derive a system of ordinary differential 
equations for the vector a(t). For this purpose, let us write ()6.7)) in the form 

Q,(Z(t), Tj(t)) = 0, j — 1, ... ,6, 0<t<U, (6.8) 

where the vectors Tj(t) = Tj(a(t)) span the tangent space Ts( a (t))S- Note that a{t) = (b(t),v(t)), 
where 

|u(*)| <u < 1, 0<t<t*, (6.9) 

by Lemma 13.41 iii). It would be convenient for us to use some other parameters (c, v) instead of 

a = (b, v), where c(t) = b(t) — \ v(r)dr and 

Jo 

c(t) = b(t) -v(t) =w{t) -v{t), 0<t<U. (6.10) 

We do not need an explicit form of the equations for (c,v) but the following statement. 

Lemma 6.2 Let Y{t) be a solution to the Cauchy problem and \4-l\j , / 1 6'. ty) hold. Then 

(c(t),v(t)) satisfies the equation 

(S) ) =A/ >W> Z W)' 0<t<t„ (6.11) 

where 

N{a,Z) = 0{\\Z\\%) (6.12) 

uniformly in o G E(ti). 

Proof We differentiate the orthogonality conditions 1)6.8)1 in t, and obtain 

= tt( y Z,T j ) + tt( y Z,i j ) = tt(AZ + T + N,T j )+tt(Z,T j ), 0<t<t*. (6.13) 

First, let us compute the principal (i.e. non- vanishing at Z = 0) term f2(T, Tj). For j = 1, 2, 3 one 
has by KTty . (B 

fi ( T ' T j) = ~ ^2(pfi(n, r i) + ^( r «+3, T j)) = ^2 fi ( r J' Tl +^ i]l = 

i i i 

where the matrix Q + is defined by ()3.6|) . Similarly, 

Q(T, r i+3 ) = - ^2(ciQ(ti, r j+3 ) + vitt(n +3 , r j+3 )) = ^ ^( r i+3, n)ci = - ^ njq. 
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n<r,T)=( ° T (v) )(i)=n { v)(i ) (6 .i4) 



As the result, we have by ([3.5)1 . 

( -O'-(r) / I r / " " li ' I r 

in the vector form. 

Second, let us compute Q(AZ, Tj). The skew-symmetry (|5.5|) implies that Q(AZ, Tj) = —Q(Z, Atj] 
Then for j = 1, 2, 3, we have by (j5.1U|) . 

Sl(AZ, tj) = -tt(Z, c ■ Vrj), (6.15) 

and similarly, 

Q(AZ, r j+3 ) = -Q(Z, c ■ Vr J+3 + Tj) = -Q(Z, c • Vr j+3 ) - Q(Z, Tj) = -Q(Z, c ■ Vr j+3 ), (6.16) 
since Q(Z, Tj) = 0. 

Finally, let us compute the last term Q(Z,fj). For j = 1, ... ,6 one has fj = b-VbTj + v -V v Tj = 
v ■ \7 v Tj since the vectors Tj do not depend on b according to ()3.3|) . Hence, 

tt(Z, fj) = tt(Z, v ■ V v Tj). (6.17) 

As the result, by ()6.14J) - (j6.17|) . the equation ()6.13|) becomes 

= Q(v) ( 6 .^j+ M (a, ^) ( • ) + M>{<r, Z), (6.18) 

where the matrix A4 (cr, Z) = (D(\\Z\\-p), and Mo{c,Z) = uniformly in a 6 £(t>) and 

11^11-/9 < r -p{v)- Then, since Q(v) is invertible by Lemma f3. 21 and H^H-^ is small, we can resolve 
(|6.18|) with respect to the derivatives and obtain equations (|6.11|) with M = (9(||Z|| 2 /3 ) uniformly 
in a G £({;). ■ 

Remark 6.3 The equations (jfi.llj) . ()6.12)1 imply that the soliton parameters cit) and v(t) are 
adiabatic invariants (see |3j). 

6.2 Decay for the Transversal Dynamics 

In Section 11 we shall show that our main Theorem 12.51 can be derived from the following time 
decay of the transversal component Z(t): 



Proposition 6.4 Let all conditions of Theorem \2. ,51 hold. Then = oo, and 

C(p,v, d ) 



-* ~ (l + \t\)^ ' ~ °- (6 - 19) 

We shall derive ()6.19|) in Sections 7-11 from our equation ()4.8|) for the transversal component Z(t). 
This equation can be specified by using Lemma f6. 21 Indeed, the lemma implies that 

\\T{t)\\p<C(v)\\Z{t)\\ 2 _ p , 0<t<U, (6.20) 
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by ()4.1Uj) since w — v — c. Thus, equation (j4.8J) becomes 



Z(t) = A(t)Z(t) + N(t), 0<t<U, (6.21) 

where A(t) = A v m >w m, and N(t) := T{t) + N(t) satisfies the estimate 

\\N{t)\\ p <C\\Z{t)\\ 2 _p, 0<t<U (6.22) 

In the remaining part of our paper we mainly analyze the basic equation (|6.21|) to establish 
the decay (|6.19|) . We are going to derive the decay using the bound ()6.22|) and the orthogonality 
condition (|6.7j) . 

Let us comment on two main difficulties in proving ()6.19|) . The difficulties are common for 
the problems studied in [5]. First, the linear part of the equation is non-autonomous, hence we 
cannot apply directly known methods of scattering theory. Similarly to the approach of j3] , we 
reduce the problem to the analysis of the frozen linear equation, 

X(t) = A 1 X(t), tEJR, (6.23) 

where Ai is the operator A Vim defined by (j4.9|) with v x = v(ti) for a fixed ti G [0,t*). Then we 
estimate the error by the method of majorants. 

Second, even for the frozen equation ()6.23|) . the decay of type (16.19)) for all solutions does not 
hold without the orthogonality condition of type (|6.7|) . Namely, by ()5.13|) equation (|6.23|) admits 
the secular solutions 

3 3 

X(t) = J2 C Mvi) +J2 D A-rj(vi)t + r j+3 ( Vl )]. (6.24) 
i i 

The solutions lie in the tangent space Ts( ai )S with a± = (pi, vi) (for an arbitrary b\ G IR) that 
suggests an unstable character of the nonlinear dynamics along the solitary manifold (cf. Remark 
14.11 ii)). Thus, the orthogonality condition (J6.7|) eliminates the secular solutions. We shall apply 
the corresponding projection to kill the unstable "longitudinal terms" in the basic equation (J6.21)) . 

Definition 6.5 i) For v G V, denote by U v the symplectic orthogonal projection of £ onto the 
tangent space Ts( a )S, and write P v = I — n„. 

ii) Denote by Z v = P V S the space symplectic orthogonal to Ts( a )S with o = (b, v) (for an arbitrary 
belR). 

Note that by the linearity, 

u v z = J2 ^d v >A v Mn(v), z), ze£, (6.25) 

with some smooth coefficients Uji(v). Hence, the projector 11^ does not depend on b (in the 
variable y = x — b), and this explains the choice of the subindex in ILj and P v . 
We have now the symplectic orthogonal decomposition 

£ = T s(fT )S + Z v , <r=(b,v), (6.26) 

and the symplectic orthogonality ()6.7|) can be represented in the following equivalent forms, 

U v(t) Z{t) = 0, P v{t) Z{t) = Z(t), 0<t<U. (6.27) 
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Remark 6.6 The tangent space Ts( a )S is invariant under the operator A v>v by Lemma 15.41 i). 
hence the space Z v is also invariant by (|5.5|) . namely: A V>V Z G Z v on a dense domain of Z G 

In Sections 12-18 below we will prove the following proposition which is one of the main 
ingredients to proving (jdl9j) . Let us consider the Cauchy problem for equation (|6.2Hj) with A = 
A VjV for a fixed v G V. Recall that the parameter /3 > 3/2 is also fixed. 

Proposition 6.7 Let M.tA) and ALIO)) hold, \v\ < v < 1, and X G S. Then 

i) The equation \6. 2ty . with Ai = A = A %v , admits a unique solution e At X := X(t) G C(JR, £) 
with the initial condition X(0) = X . 

ii) For X G Z v n S/3, the solution X(t) has the following decay, 

\\e At X \\-p< (il^ /2 \\ X °h> te]R - ( 6 - 28 ) 

Remark 6.8 The decay is provided by two fundamental facts which we will establish below: 

i) the null root space of the generator A coincides with the tangent space Ts( a )S, where a = (b, v) 
(for an arbitrary b G IR), and 

ii) the spectrum of A in the space Z v is purely continuous. 

7 Frozen Form of Transversal Dynamics 

Now let us fix an arbitrary t\ G [(),£*), and rewrite the equation (|6.21j) in a "frozen form" 

Z(t) =A 1 Z(t) + (A(t)-A 1 )Z(t) + N(t), 0<t<t m , (7.1) 
where A 1 = A v{tl)Atl) and 



A(t) - A, 



( [w(t)-v(ti)] ■ V \ 

[w(t)-v(tx)] ■ V 

B v (fi — B v ( tl -) 

V o o (V(Vyt)-Vy tl )),V/>) o J 



The next trick is important since it enables us to kill the "bad terms" [w(t) — v(ti)] ■ V in the 
operator A(t) — A\. 

Definition 7.1 Let us change the variables (y,t) i— ► (yx,t) = (y + dx(t),t) where 

d x {t) := / (w(s) -v(ti))ds, < t < t x . (7.2) 

Next, let us write 



Z 1 {t) = {*{ yi - dt{t), t), U( yi - d 1 (t),t), Q(t), P(t)). (7.3) 
Then we obtain the final form of the "frozen equation" for the transversal dynamics 

Z x {t) = A x Z x (t) + B 1 (t)Z 1 (t) + Ni(t), < t < t 1: (7.4) 
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where Ni(t) = N(t) expressed in terms of y — y\ — dx(t), and 
Bi(t) 



( \ 



B v (t) — B v (t x ) 

\ (V(^ w -^ (<l) ),Vp) / 



At the end of this section, we will derive appropriate bounds for the "remainder terms" Bi(t)Zi(t) 
and Ni(t) in (|7.4|) . First, note that we have by Lemma fo.2\ 

\B v[t) - B v{tl) \ < | [ v(s) ■ V v B v{s) ds\ < C [ \\Z{s)\\%ds. (7.5) 

Jtl Jt 

Similarly, 

|(V(^ ( t) ~^(*i)),Vp>| <C ^ \\Z(s)\\ 2 _ p ds. (7.6) 
Let us recall the following well-known inequality: for any a G IR 

(l + \y + x\T < (1 + \y\) a (l + x, y G IR 3 . (7.7) 



Lemma 7.2 For (\l/, IT, Q, P) G £ Q icit/i any a G IR i/ie following estimate holds: 

\\(*(y 1 -d 1 ),Il(y 1 -d 1 ),Q,P)\\ a < ||(*,n,Q,P)|| a (l + |di|) |a| , ^ G IR 3 . (7.8) 
Proof Let us check the estimate only for one component, say, for IT. One has by (J7.7|) 

l|n(jfc-di,t)||iU = / \Il(y 1 -d 1 ,t)\ 2 (l + \y 1 \) 2a dy 1 = J \U(y, t)| 2 (l + \y + d 1 \) 2a dy < 

\U(y,t)\ 2 (l + \y\) 2a (l + Id^dy < (1 + {d^ \\U\\l a , 
and the lemma is proved. ■ 
Corollary 7.3 The following bound holds 

\\Ni{t)\\ p < (1 + {dMYmW-p > < t < t x . (7.9) 
Indeed, applying the previous lemma, we obtain from ()6.22|) that 

\\N X {t)h < (i + \d x m p \m,z(t))h < (i + m^mW-p- 

Corollary 7.4 The following bound holds 

PM*)2i(*)ll/J < C\\Z(t)\\- P J* \\Z{r)\\%dT , < t < h. (7.10) 

For the proof we apply Lemma 17.21 to ()7.5|) and ()7.6|) and use the fact that Bi(t)Zi(t) depends 
only on the finite-dimensional components of Z\{t). 
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8 Integral Inequality 



Equation (|7.4jl can be represented in the integral form: 

Z x (t) = e Al %(0) + [ e Al(t-a) [SiZi(s) + Ni(s)]ds, < t < t x . (8.1) 
Jo 

We apply the symplectic orthogonal projection Pi := P v tti) to both sides, and get 

PxZ^t) = e^P^O) + [ e Al( t-''*P 1 [B 1 Z 1 (s) + N^sflds. 

Jo 

We have used here that P\ commutes with the group e Alt since the space Z\ := P\S is invariant 
with respect to e Alt by Remark 16.61 Applying ()6.28|) we obtain 



\PlZ! 



(1 + t) 3 / 2 ' 



o (l + \t-s\)^ 



\PtlBxZtis) + NxisMpds. (8.2) 



The operator Pi = I — Hi is continuous in £p by ()6.25|) . Hence, from (|8.2|) and ()7.9|) . ()7.10|h we 
obtain 



l^i(*)ll-/3<7^£ll^(0)H. 



+C(di(t)) 



(1 + t) 3 / 2 
* 1 



|Z( S )|L 



Z(r)|| 2 _^+ 11^)11% 



'o (l + |t-3|) 3 / 2 
where d\(t) := sup 0<s<< |g?i(s)|. 
Definition 8.1 Let be the exit time 

t[ = sup{t G [0,t*) : di(s) < 1, < s < t}. 

Now (fOjt implies that for t x < ^ 

C 



ds, 0<t< t u (8.3) 



5.4) 



Vo (l + \t-s\)^ 



\Z(s) 



\Z(r)\t R dr+\\Z(s)\\ 2 _ 



ds, 0<t< ti. 



J.5) 



9 Symplectic Orthogonality 



Finally, we are going to change P\Z\{£) by Z(t) in the left hand side of ()8.5|) . We shall prove that 



this change is possible indeed by using again the smalness condition ()2.12|) . For the justification 
we reduce the exit time further. First, introduce the "majorant" 



m(t) := sup(l + S ) 3 / 2 ||Z( S )||_ /3 , 
se[o,t] 



t E [0,U 

Denote by e a fixed positive number (which will be specified below). 



(9.1) 
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Definition 9.1 Let t'l be the exit time 

t'l = sup{t G [0, t'J : m{s) <e, < s < t}. (9.2) 

The following important bound ()9.3j) enables us to change the norm of PiZ\(t) on the left 
hand side of (|8.5|) by the norm of Z(t). 

Lemma 9.2 For sufficiently small e > 0, we have 

\\Z(t)\\-f3 < C\\PiZ 1 {t)\\- Pl 0<t<h, (9.3) 

for any ti < t", where C depends only on p and v. 

Proof The proof is based on the symplectic orthogonality ()6.27j) . i.e. 

U v(t) Z(t) = 0, te[0,h}, (9.4) 

and on the fact that all the spaces Z{t) := P v (t)£ are almost parallel for all t (see Fig. 2). Namely, 





Figure 2: Symplectic orthogonality. 

we first note that < C||^iWII-/3 by Lemma 17^1 since |di(t)| < 1 for t < h < t'l < t' # . 

Therefore, it suffices to prove that 



\\Z 1 (t)\\-p<2\\PiZ 1 
This estimate will follow from 

1 



< t < t x . 



Wn^z^t)^ < -wz^t)^, o<t<t u 

since P\Zi{t) = Z\{t) — ILj^Zi^i). To prove (|9.fjjl . we write (|9.4jl as 

U^Z^t) = 0, te[0,h], 



(9.5) 



(9.6) 



(9.7) 



where Tl^^Zxit) is Yl v ^Z{t) expressed in terms of the variable yi = y + d\(t). Hence, (|9.6|) 
follows from ()9.7|) if the difference LLj^) — LLj^y is small uniformly in t, i.e. 



|n„( tl ) — LLj^xH < 1/2, 



<t <t x . 



(9. 
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It remains to justify ()9.8|) for any sufficiently small e > 0. We will need the formula ()6.25|) and 
the following relation which follows from ()6.25|) : 

U^xZtit) = ^n i ,(u(t))r i ,i(i;(t))n(7i, 1 (T;(0),^ 1 (t)), (9.9) 
where Tj t i{v{t)) are the vectors Tj(v(t)) expressed via the variables y\. In detail (cf. (|3.3|) ). 

j = 1,2,3, (9.10) 



T jA v ) '■= (-9jtpv(yi - ^(t)), -d j ir v (y 1 - di(t)),ej,0), 
Tj+3,i(v) := {d Vj i) v (yi - d v .TT v (yi - di(t)), 0, d Vj p v ), 



where v = v(t). Since < 1, and the functions Vr.,- are smooth and rapidly decaying at 

infinity, Lemma f7.2l implies that 

IMu(t)) - r^*))^ < Cl^Ct)!^, < t < ti (9.11) 

for all j = 1, 2, . . . , 6. Furthermore, 

Tj(v(t)) - Tj(v(ti)) = v(s) ■ V v Tj(v(s))ds, 

and therefore t 

\\rMt))-rMti))h<C J t ' \v{s)\ds, < t < fc. (9.12) 

Similarly, 

in,-, («(*)) - n^ti))] = | / 1 v( s ) ■ v„n i ,(«(s))ds| < c [ 1 |t>(s)|ds, o < t < ti, (9.13) 

since \V v Ilji(v(s)) \ is uniformly bounded by (|f).9)l . Hence, the bounds ()9.8|1 will follow from (|fi.25|) . 
(|9.9|) and ()9.11|) - (|9.13|l if we shall prove that \di(t) \ and the integral on the right hand side of (J9.12|) 
can be made as small as desired by choosing a sufficiently small e > 0. 
To estimate di(t), note that 

w(s) - v(ti) = w(s) - v(s) + v(s) - v(ti) = c(s) + f v(r)dT (9.14) 

J s 

by (|6.10|) . Hence, the definitions (|7.2|) . ()9.1|) . and Lemma f6.2l imply that 
|di(t)| = | f (w(s) - < [ 1 (\c(s)\ + [ 1 \v{r)\dr) ds 

Jt! Jt V </s / 

< Cm 2 ^) jf 1 (r^T^ + / 1 (TT^) dS ~ ° m2{tl) ~ C£2 ' ° - t - (9 ' 15) 
since t\ < t". Similarly, 

/ \v(s)\ds ^Cm 2 ^) b —<Ce 2 , 0<t<h. (9.16) 

Jt Jt (i + sr 

The proof is completed. ■ 



24 



10 Decay of Transversal Component 



Here we prove Proposition 16.41 

Step i) We fix e > and t" = t"(e) for which Lemma E21 holds. Then a bound of type (|8.5jl holds 
with ||P 1 Z 1 (t)||_^ replaced by on the left hand side: 



< 



C 



+c 



(1 + tf/ 2 
1 



1^(0)11/ 



\Z(s\\\- t 



Z{r)\\%dr+\\Z{s)f_ 



ds, < t < h 



'Q (l + |i- S |) 3 / 2 

for t x < ^. This implies an integral inequality for the majorant 

m(t) := sup {l + sf' 2 \\Z{s)\\^ . 
se[o,t] 

Namely, multiplying (jlO.lj) by (1 + t) 3 ^ 2 and taking the supremum in t G [0, ti], we get 

m(s) 



(10.1) 



m(t 1 )<C||Z(0)|| /3 + C' sup / ^ T ) h3/2 

te[o,*i] Jo U + l r _ s \) ' 



h m 2 (r)dT m 2 (s) 



_(l + s )3/2 J g (l + r ) 3 ' (1 + s^ 3 
for t\ < t". Taking into account that m(t) is a monotone increasing function, we get 



ds 



m{h) < C\\Z{G)\\ P + Cim^h) + m 2 (t!)]/(t 



h < c 



10.2) 



where 



(1 + t) 3 / 2 

t!Z]Jo (i + \t-s\y/ 2 



dr 



+ 



_(l + s) 3 / 2 J s (l + r) 3 ' (1 + s) 3 J 
Therefore, ()10.2j) becomes 

m(h) <C||Z(0)||^ + C'7[m 3 (t 1 )+m 2 (t 1 )], ^ < C 
This inequality implies that m(t\) is bounded for t\ < t", and moreover, 



ds < I < oo, £i > 0. 



u < t: 



(10.3) 



(10.4) 



since m(0) = ||Z(0)||/3 is sufficiently small by (|H.lljl . 

Step ii) The constant C*i in the estimate ()10.4|) does not depend on t*, ^ and t" by Lemma E21 
We choose a small d in (l2~T2l such that 11^(0)1^ < e/{2C x ). This is possible by (I3~TT1) . Then the 
estimate ()10.4|) implies that £" = ^ , and therefore ()10.4|) holds for all t\ < t'^. Then the bound 
f!9.15|) holds for all t < t*. Choose a small £ such that the right hand side in ()9.15|) does not exceed 
one. Then t'^ = t*. Therefore, (|10.4|) holds for any t\ < t" = t*, hence (|6.4|) also holds if H-Z^O)^ is 
sufficiently small. Finally, this implies that = oo. Hence we also have t" = ^ = 00, and (|10.4j) 
holds for any t\ > if d is sufficiently small. ■ 
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11 Soliton Asymptotics 

Here we prove our main Theorem 12.51 under the assumption that the decay ()6.19j) holds. Let us 
first prove the asymptotics ()2.13j) for the vector components, and then the asymptotics (j2.14j) for 
the fields. 

Asymptotics for the vector components It follows from ()4.3|) that q = b + Q, and from (j6.21j) . 
fdb . and (HHJ) that Q = B v(t) P + 0(\\Z\\ 2 _ p ). Thus, 

q = b + Q = v(t) + c(t) + B v[t) P(t) + 0(\\Z\\ 2 _ p ). (11.1) 
Equation (J6.11|) together with estimates (|6.12j) and (J6.19j) imply that 

|c(t)| + |»(*)| < C -^0 1 , t > 0. (11.2) 

Therefore, c{t) = c + + 0{t~ 2 ) and v (i) = v + + 0{t~ 2 ), t — > oo. Since \P\ < \\Z\\^p, the estimate 
()6.19|) together with relations (|11.2|) and (jll.lj) imply that 

q{t)=v + + 0{r" 2 ). (11.3) 

Similarly, 

b(t)=c(t)+ v(s)ds = v + t + a + + 0(t- 1 ), (11.4) 
Jo 

and hence the second part of ()2.13|) follows: 

q(t) = b(t) + Q(t) = v + t + a + + Oir 1 ), (11.5) 
since Q(t) = £>(r 3 / 2 ) by (foTHjl . 

Asymptotics for the fields We apply the approach developed in fSJ|23]- ^ or ^ ne P ar ^ °^ 
the solution, F(t) = (ip(x,t),TT(x,t)), in the original variable x, let us define the accompanying 
soliton field as F v ^(t) = (^ v (t)(a ; — q(t)), ^v(t)( x — where we now set v(t) = q(t), cf. (|1 1. 1|) . 

Then for the difference Z(t) = F(t) — F y M(t) we obtain easily, from the first two equations of the 
system (jl.2j) . the inhomogeneous Klein-Gordon equation (221 (2.5)], 

Z(t) = A Z(t)-v-V v F v(t) (t), A ^,n) = (7r,(A-m 2 )^). 

Then 

Z(t) = W o (t)Z(0) - [ W (t - s)[v(s) • V v F v(s) (s)}ds, (11.6) 
Jo 

where Wo(t) is the dynamical group of free Klein-Gordon equation. To obtain the asymptotics 
(|2~T4l it suffices to prove that Z(t) = W (t)^f + + r + (t) for some * + 6f and that ||r + (t)||^ = 
0{ir x l 2 \ This is equivalent to the asymptotics 

W (-t)Z(t) = * + + r' + (t), \\r' + (t)\\f = Oir 1 / 2 ), (11.7) 

since Wo(t) is a unitary group on the Sobolev space JF by the energy conservation for the free 
Klein-Gordon equation. Finally, the asymptotics ()11.7|) hold since (I11.6J) implies that 

W (-t)Z(t) = Z(0) - f W (-s)R(s)ds, R(s) = v(s)-V v F v(s) (s), (11.8) 

Jo 
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where the integral on the right hand side of (jll.8|) converges in the Hilbert space T with the 
rate 0{t~ 1 / 2 ). The latter holds since || W (-s)R(s) \\ F = 0(s~ 3/2 ) by the unitarity of W (-s) and 
the decay rate ||i?(s)||^ = 0(s~ 3 ! 2 ). Let us prove this rate of decay. It suffices to prove that 
|v(s)| = C(s~ 3/2 ), or equivalent^ \p(s)\ = 0(s~ 3/2 ). Substitute (JOD to the last equation of (TT2j) 
and obtain 



P(t) = I fy v ( t )(x-b(t)) + V(x-b(t),t)]Vp(x-b(t)-Q(t))dx (11.9) 
^ {t) {y)V p{y)dy + [ 4> v[t) {y)[Vp{y-Q{t))-Vp{y)]dy+ [ *{y,t)Vp{y - Q{t))dy. 



The first integral on the right hand side is zero by the stationary equations (J2.7j) . The second 
integral is 0(t~ 3 ^ 2 ), which follows from the conditions f)1.9|) on p and the asymptotics Q(t) = 
0(t~ 3 / 2 ). Finally, the third integral is (D(t~ 3 / 2 ) by estimate (|6.19|) . This completes the proof. ■ 



12 Decay for the Linearized Dynamics 

In remaining section, we prove Proposition 16 . 71 to complete the proof of the main result (Theorem 
12. 5 j) . Here we discuss the general strategy of proving the Proposition. We apply the Fourier- 
Laplace transform 

POO 

X(X) = e- xt X(t)dt, Re A > (12.1) 
Jo 

to (|fi.23)l . According to Proposition 16.7) we can expect that the solution X(t) is bounded in the 
norm || • ||_^. Then the integral (j!2.1j) converges and is analytic for ReA > 0, and 

\\X(X)\\^<^, ReA>0. (12.2) 

Let us derive an equation for X(X) which is equivalent to the Cauchy problem for (jd23j) with the 
initial condition X(0) = X G £_p. We shall write A and v instead of Ax and vi in all remaining 
part of the paper. Applying the Fourier-Laplace transform to ([6.23)1 . we get that 

XX(X) = AX(X) + X , Re A > 0. (12.3) 

Let us stress that ()12.3|) is equivalent to the Cauchy problem for the functions X(t) G C&([0, oo); £_^) 
Hence the solution X(t) is given by 

X(X) = -(A- A) _1 X , Re A > (12.4) 

if the resolvent R(X) = (A — A)^ 1 exists for Re A > 0. 

Let us comment on our following strategy in proving the decay ()6.19|) . We shall first construct 
the resolvent R(X) for ReA > and prove that this resolvent is a continuous operator on £_p. 
Then X(X) G S-p and is an analytic function for ReA > 0. After this we must to justify that 
there exists a (unique) function X(t) G C([0, oo); £-p) satisfying ()12.1j) . 

The analyticity of X(X) and the Paley- Wiener arguments (see |22j) should provide the existence 
of a £-p - valued distribution X(t), t G 1R, with a support in [0, oo). Formally, 

X(t) = — [ e itJt X(iu + 0)du, t G H. (12.5) 
2tt J n 
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However, to establish the continuity of X(t) for t > 0, we need an additional bound for X(iu + 0) 
for large values of |a>|. Finally, for the time decay of X(t), we need an additional information on 
the smoothness and decay of X(iu + 0). More precisely, we must prove that the function X(iu + 0) 
has the following properties: 

i) it is smooth outside u = and u = ±/i, where \i = /z(t>) > 0, 

ii) it decays in a certain sense as \u\ — > oo, 

hi) it admits the Puiseux expansion at u = ±p, 

iv) it is analytic at u> = if X G Z v := P V S and X G £p. 

Then the decay (|6.19|) would follow from the Fourier-Laplace representation (|12.5|) . 

We shall check with detail properties of the type i)-iv) only for the last two components Q(X) 
and P(A) of the vector AT (A) = (tf (A), 6(A), Q(X), P(A)). The properties provide the decay 
for the vector components Q(t) and P(t) of the solution X(t). 

However, we will not prove the properties of the type i)-iv) for the field components A) 
and n(a;, A). The decay (j6.19|) for the field components is deduced in Section 18 directly from the 
time-dependent field equations of the system (J6.23|) . using the decay of the component Q(t) and 
a version of strong Huygens principle for the Klein-Gordon equation. 



13 Constructing the Resolvent 

To justify the representation (jl2.4J) . we construct the resolvent as a bounded operator in £_p for 
Re A > 0. We shall write U(y), Q, P) instead of (y, A), fl{y, A), Q{\), P(A)) to simplify 

the notations. Then (|12.3j) reads 



(A-X) 



n 

Q 
\P ) 



( *o \ 
n 

Qo 



where A 



( * \ 

n 

Q 



( n + v ■ w \ 

A\l> - m 2 ^ + v ■ Vn + Q ■ Vp 

B V P 

\ -(v^,p) + (v^,g-v P ) j 



This gives the system of equations 

Tl(y) + v ■ V*(y) - W(y) = ~Mv) 

AV(y) - m 2 ^{y) + v ■ VU{y) + Q ■ V 'p(y) - An(y) = -n (y) 

S„P - AQ = -Qo 

-<V*(y),p(y)) + (VMv),Q- Vp(y)> - AP = -P 
S'tep ^ Let us study the first two equations. In the Fourier space they become 

n(jfe) - ivjfe*(jfe) - A*(jfe) = -*oO) 

(-A; 2 - m 2 )^(k) - (ivk + A)fl(A;) = — n (A;) + iQkp(k) 
Let us invert the matrix of the system and obtain 

[(ivk + A) 2 + k 2 + m 2 ]- 1 



yeJR 3 



;i3.1) 



:i3.2i 



-(ivA; + A) 1 
-{k 2 +m 2 ) -(ivk + X) 



-i 



-(ivfc + A) 
k 2 + m 2 



-1 

(ivk + A) 



28 



Taking the inverse Fourier transform, we obtain the corresponding fundamental solution 

Gx(y)=( % V ~J " X ,)gx(v), (13.3) 



-A + m 2 v • V - A f 

where g\(y) is the unique tempered fundamental solution of the determinant 

D = D(X) = -A + m 2 + (-v ■ V + A) 2 . (13.4) 

From now on we use the system of coordinates in x-space in which v = (\v\, 0, 0), hence vk = \v\ki, 
and 

1 1 

9x(y) = F k \ k2 + m2 + ^ - A)2 = F-_+ y - + m2 + - A)2 , y E IR 3 . (13.5) 

Note that the denominator does not vanish for Re A > 0. This implies 

Lemma 13.1 The operator G\ with the integral kernel G\{y — y') is continuous as an operator 
from H\]R 3 ) © L 2 (R 3 ) to H 2 (1R 3 ) © H^R 3 ) for Re A > 0. 

Thus, formulas ()13.2|) and ()13.3|) imply the convolution representation 
$ = -{y. V-A)^ A **o + ^A*n + (^A* Vp)-Q 

IT = -(-A + m 2 )g x *^ -(vV - \)g x *Tl -(v V - \){g x * Vp) • Q 



(13.6) 



Step ii) Let us compute g\(y) explicitly. First consider the case v — 0. The fundamental solution 
of the operator —A + m 2 + A 2 is 

e ~ K \y\ 

gx(y) = rTT , (13.7) 

where 

K 2 = m 2 + A 2 , Re « > for Re A > 0. (13.8) 
Thus, in the case t> = we have 



Gx(y-y') 



-A + m 2 -A i 47r|?/-y / | 



For general v = (\v\,0,0) with \v\ < 1 the denominator in ()13.5|) . which is the Fourier symbol of 
D, reads 

D(k) = k 2 + m 2 + (i\v\ki + A) 2 

= (1 - v 2 )k\ + k 2 2 + k 2 z + 2i\v\h\ + \ 2 + m 2 

= (l-v 2 )(k l + ^-) 2 + k 2 2 + kl + K 2 , (13.9) 
1 — tr 



where 



K 2 = + a 2 + m 2 = + m 2 . (13.10) 

1 — v A 1 — v z 
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Therefore, setting 7 := l/vl — v 2 , we have 

k = 7a/ A 2 + p 2 , p := m/y. (13.11) 

Return to x-space: 

D = — ^r(Vi + 7^i) 2 — V 2 . — V 2 + K 2 , Ki:=7HA (13.12) 

7 

Define y\ := yy\ and Vi := d/dyi. Then 

D = _(Vi + ki) 2 - V 2 . - V 2 + k 2 . (13.13) 
Thus, its fundamental solution is 



e 



-K\V\-KlVl 



9x(y) = — j-pj — > 1/ := 2/2, 2/3), (13.14) 

47T|j/| 

where we choose Re k > for Re A > 0. Let us note that 

< Re K\ < Re k, Re A > 0. (13.15) 

This inequality follows from the fact that the fundamental solution decays exponentially by the 
Paley- Wiener arguments since the quadratic form ()13.9j) does not vanish in a complex neigh- 
borhood of the real space IR 3 for Re A > 0. Let us state the result which we have got above. 



Lemma 13.2 i) The operator D = D(X) is invertible in L 2 (IR 3 ) for Re A > and its fundamental 
solution \1S. 14]) decays exponentially. 

ii) Formulas \1S. 14]) and H13.ll)) imply that, for every fixed y, the Green function g\(y) admits 
an analytic continuation (in the variable X) to the Riemann surface of the algebraic function 
■\/A 2 + p 2 with the branching points ±ifi. 

Step Hi) Let us now proceed with the last two equations ()13.1|) . 

-XQ + B V P = -Q , (VW Q ■ Vp) - (V*, p) - XP = -P . (13.16) 

Let us eliminate the field \1/ by the first equation (|13.6|) . Namely, rewrite the equation in the form 
= + ^ 2 (^o, n ), where 

MQ) = Q-(gx*Vp), tf 2 (*o,n o ) = -(u-v-A)0A**o + 0A*no. (13.17) 

Then we have 

(v* jP ) = (w x (Q),p) + (v* 2 (* ,n ),p), 

and the last equation in (|13.16|) becomes 

(Wv, q ■ Vp) - (VMQ),p) - ap = -p + <w 2 (* , no), p) =■ -K 
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Let us first compute the term (Vip v , Q ■ Vp) = J^ji^^v, Qjdjp) = ^2j(ViJ) v , djp)Qj. Applying the 
Fourier transform F y ^i the Parseval identity, and (jA.5|) we see that 



(diip v ,djp) = (-iki$ v (k),-ikjp(k)) = (ki$ v (k),kjp(k)) 



;,*,«*)> = " / J^f^t, = -K V . (13.18) 



"k 2 + m 2 -{\v\kx) 2 ' Jrv y/ 7 k 2 + m 2 - (\v\hr) 

As the result, (Vip v ,Q ■ Vp) = —KQ, where K is the 3x3 matrix with the matrix elements 
Kij. The matrix K is diagonal and positive definite since p(k) is spherically symmetric and not 
identically zero by (|1.1U|) . 

Let us now compute the term — (V\I/i,p) = (^i, Vp). We have 

d t p) = C^2(g\ * djp)Qj, d^) = ^2(g x * djp, dip)Qj = ^ %(A)Q j 

J 3 3 

since = Q • (g\ * Vp), and by the Parseval identity again, we have 
H i:j (X) : = (g x * djp, dip) = (igx{k)kjp(k),ikip(k)) 

%km .»M(*0) = f -r, M/KWfk (13ig) 



k 2 + m 2 + (i\v\ki + X) 2 ' J k 2 + m 2 + (i\v\ki + X) 2 

The matrix H is well defined for Re A > since the denominator does not vanish (or g\(x) 
exponentially decays). The matrix H is diagonal similarly to K. Indeed, if i ^ j, then at least 
one of these indexes is not equal to one, and the integrand in ()13.19|) is odd with respect to the 
corresponding variable. 

As the result, — (V\l/i,p) = HQ, where H is the diagonal matrix with matrix elements Hjj, 
1 < j < 3. Finally, the equations (jl3.16j) become 

M(A) (?) = (??). ""ere M(A) = ( R \ x) ) , (13.20) 

where the matrices K and H(X) are diagonal. 

Step iv) Assume for a moment that the matrix M(X) is invertible for Re A > (later we shall 
prove that this the case indeed). Then 

f^)=M-W Re A > 0. (13.21) 

Finally, formulas (|13.21|) and (|13.6j) give the expression of the resolvent R(X) = (A-Xy 1 , Re A > 0. 

Lemma 13.3 The matrix function M(A) (M~ 1 (\)) admits an analytic (meromorphic) continua- 
tion from the domain Re A > to the Riemann surface of the function a/A 2 + p? . 

Proof The analytic continuation of M(A) exists by Lemma [13. 21 ii) and the convolution expres- 
sions in (|13.19|) since the function p(x) is compactly supported by (|1.9|) . The inverse matrix is then 
meromorphic since it exists for large Re A: this follows from ()13.2U|) since H(X) —>■ as Re A — > oo 

by (gni - ■ 
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14 Analyticity in the Half-Plane 



Here we prove the following proposition. 

Proposition 14.1 The operator-valued function R(X) : £ 

Proof It suffices to prove that the operator A — X : £ - 
Re A > 0. Recall that A = A vv where \v\ < 1. 



> £ is analytic for Re A > 0. 
£ has bounded inverse operator for 



Step i) Let us prove that Ker (A — A) =0 for Re A > 0. Indeed, assume that the vector X\ = 
Tlx, Q\i P\) £ £ satisfies the equation (A — X)X\ = 0, that is X x is a solution to ()13.1j) with 
\I>o = n = and Q — P — 0. We must prove that X x = 0. 

Let us first show that P\ = 0. Indeed, the trajectory X := X x e xt G C(JR,£) is the solution to 
the equation X = AX of type (|5.1|) with w = v. Then H V)V (X(t)) grows exponentially by (|5.14j) . 
since the matrix B v is positive. This growth contradicts to the conservation of TC V , V , which follows 
from Lemma 15.11 ii) because X(t) G C 1 (IR, £ + ). The latter inclusion follows from Lemma Il3.ll 
since Tlx) satisfies equations p3.(ij) with ^ = n = and Q = Qx- 

We now have XQx = B v Px = by the third equation of ()13.1|) . and hence Qx = because 
A ^ 0. Finally, ^ A = 0, n A = by equations (flTH|l with Q = Q x = 0. 

Step ii) One has 



(A-X) 



( * ^ 




n 




Q 




\p J 





v-V-X 
A — m 2 




1 

V 





-A 









\ 




•Vp 







n 


-A 


By 




Q 


b vi - Vp) 


-X 




\p J 



Thus, A — X 



An 



X 


= ^4o + T, 


where 








I 


v-V-X 


1 








\ 




A — m 2 


v-V-X 


















-A 







\ 











-A 


/ 



T 

















\ 










■Vp 


















B v 




{(■ 


,Vp) 





(V^, - Vp) 





/ 



The operator T is finite-dimensional, and the operator A 1 is bounded on £ by Lemma 113.11 
Finally, A - X = A Q (I + A 1 T), where A^T is a compact operator. Since we know that Ker (I + 
Aq 1 T) = 0, the operator (/ + Aq 1 T) is invertible by the Fredholm theory. ■ 



Corollary 14.2 The matrix M(X) of hi H. 2(^1 is invertible for Re A > 0. 



15 Regularity on the Imaginary Axis 

Next step should be an investigation of the limit values of the resolvent R(X) at the imaginary 
axis A = iu, u G IR, that is necessary for proving the decay (|6.19|) of the solution X(t) = 

(*(t),n(t),Q(t),p(t)). 

Let us first describe the continuous spectrum of the operator A = A VjV on the imaginary axis. 
By definition, the continuous spectrum corresponds to uj G IR such that the resolvent R(iu + 0) 
is not a bounded operator on £. By the formulas (jlH.fij) . this is the case if the Green function 
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9\{y ~ y') foils to have the exponential decay. This is equivalent to the condition that Reft = 0, 
where k is given by (jl3.11|) : k = jy/ /jp — u 2 . Thus, iu belongs to the continuous spectrum if (cf. 

(HSHD) 



\uj\> [i = myl — 



By Lemma [13.3^ the limit matrix 



M H :=M(, + 0)=(,^ +() we H, (15.1) 

exists, and its entries are continuous functions of u G IR, smooth for |u;| < /x and |u;| > /i. Recall 
that the point A = belongs to the discrete spectrum of the operator A by Lemma 15.41 i). and 
hence Miiio + 0) is (probably) not invertible either at uj = 0. 

Proposition 15.1 Let lil.ty) and M.1U\) hold, and \v\ < 1. Then the limit matrix M(iu + 0) is 
invertible for w ^ 0, w 6 E. 

Proof Let us consider the three possible cases < \uj\ < ji, \uj\ = //, and \cu\ > \l separately. Let 
us remind that the matrices K and H are diagonal with the entries 



k 2 \p{k)\ 2 dk 
k 2 + m 2 — (\v\ki 



/■ k 2 \p(k)\ 2 dk 

h ^ )= JwtJtw%tw' ReA>0 - (15 - 3) 

and i^22 = #33- Since v = (\v\, 0, 0), the matrix is also diagonal: 

/ z/ 3 \ 

B v := v{E-v®v) = i/ (15.4) 




since z/ 2 := 1 — v 2 . Let us denote -F(^) := —K + i7(iu; + 0) which is also diagonal, and let 
F\\ := F n (u), and F± := F 22 (uj) = F 33 (u). Then by (fT5~T| 

det Af(iw) = det f 7^ J = -(cu 2 + z/ 3 ^,)^ 2 + z/F ± ) 2 , u; G IR. (15.5) 

The formula for the determinant is obvious since both matrices F(uS) and B v are diagonal, hence 
the matrix M(iu) is equivalent to three independent matrices 2x2. Namely, let us transpose the 
columns and rows of the matrix M(iu) in the order (142536). Then we get the matrix with three 
2x2 blocks on the main diagonal. Therefore, the determinant of M(iu) is simply product of the 
determinants of the three matrices. 

I. First, let us consider the case < \u\ < p. Then the invertibility of M{iuj) follows from ()15.5j) 
by the following lemma. 

Lemma 15.2 For < |u;| < fi, the matrix F(u) is positive definite, i.e. Fjj(u) > 0, j = 1,2,3. 
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Proof First, let us check that the denominator in ()15.3|) is positive for A = iu with \u\ < p. 
Indeed, it equals m 2 + k 2 — (u+ H&i) 2 and we have to prove that m 2 + k 2 > oj 2 + 2u\v\ki + v 2 k\. By 
the condition \u\ < p, — myl — v 2 , it suffices to prove that m 2 + k 2 > m 2 (l — v 2 ) + 2uj\v \k\ + v 2 k\ . 
This is equivalent to k 2 + k 2 + m 2 v 2 + kf(l — v 2 ) > 2m\v\ki\/l — v 2 , which is evidently true. Thus, 

=Jmm 2 dk ( m2+k2 _\ lvlki+u)2 - m2+fc2 ^ (Hfcl)2 ) > 3 = 1,2,3. 

Let us prove that Fjj(u>) > 0. Indeed, since p(k) = p(—k), we obtain that 

2i ol I I 



Fjj (u) = dk 2 dk 3 k 2 \p(k)\ 2 ( -T-TT, TUu-r^i + 



m 2 + k 2 — (\v\k\ + lo) 2 m 2 + k 2 — (\v\ki — to) 2 

2 



m 2 + k 2 — (|f Ifcx) 2 



dfci. (15.6) 



Now it suffices to prove that the expression in brackets is positive (or positive infinite) under the 
conditions 

\v | < 1, <\lo\ < p — my/ 1 - v 2 . (15.7) 

This is proved in Appendix B. ■ 
II. uj = ±p. For example consider the case u = p. Then formula (|15.3|) reads (see (|13.9|0 : 

k 2 \p(k)\ 2 dk 



Hjjiip) 



k% + k 2 + {yki — m|v|) 5 



Now the integrand has a unique singular point. The singularity is integrable, and therefore the 
terms Fjj(p) are finite. Furthermore, the terms are positive by the integral representation (|15.6J) 
again. Hence, the matrix M(ip) is invertible. 

III. \u\ > p. Here we apply another arguments: the invertibility of M{iuj) follows from (J15.5)) by 
the methods used in j50J Chapter VII, formula (58)]. 

Lemma 15.3 If U.1U\) holds and if \u\ > p, then the imaginary part of the matrix - — :F{uS) is 



\0J\ 



negative definite, i.e. - — :lmFjj(uj) < 0, j = 1,2,3. 

\uj\ 

Proof Since F(u) = —K + H(iu + 0) where the matrix K is real, it suffices to study the matrix 
H{iu + 0). For e > 0, we have 

f k 2 \p(k)\ 2 dk 

H jj (iu + e)= / -= - 2 , 2 3 , ,/' — -, J = 1,2,3. (15.8) 

" y ' J k\ + k\ + k\ - (\v\fa + cu - ie) 2 + m 2 v ' 

Consider the denominator 

D £ (k) = k 2 + m 2 - (\v\k! +u- ie) 2 . 

It was shown above that Do{k) ^ if \u\ < p, and Do(k) vanishes at one point if \u\ = p. On the 
other hand, for \u\ > p the denominator Do(k) vanishes on the ellipsoid 

T^ik-.iuh-^ + kl + k^R 2 :=^^}, 
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where v = yl — v 2 . We shall show below that the Plemelj formula for C 1 -functions implies that 

u f k 2 \p(k)\ 2 

JmH ii (iu + 0) = — ! -rn 1 dS, (15.9) 



\U!\ 



where dS is the element of the surface area. This immediately implies the statement of the Lemma 
since the integrand in ()15.9j) is positive by the Wiener condition (jl.K)j) . 

Let us justify ()15.9|) for uj > p > (the case uj < —p < can be treated similarly). Let 
C e Q°(IR 3 ) be a nonnegative cut off function equal one when \Do(k)\ < 5 and vanishing when 
|-Do(^)| > 25. We fix a small 5 and split the integral f)15.8j) in two parts: with the factor £ in the 
integrand and with the factor 1 — (. The limit of the second term as e — > is real. Hence, we 
have to calculate the imaginary part only for 

(/r\ f k 2 \p(k)\ 2 dk 

Hfhiu + 0) = lim / ((k) M . (15.10) 

n ^oj D £ (k) 



Denote a(k) = \/k 2 + m 2 and b(k) = \v\ki + uj. Then 

1 1 1 



+ — — . (15.11) 



D £ {k) a 2 — {b — is) 2 2a(a — b + ie) 2a(a + b — is 

Note that D Q {k) ^ if b{k) = 0. Thus, b{k) ^ on T w , and therefore b{k) ^ on the support of 
C if 5 <C 1. Since b(k) > when v — {y — 1), we get that b(k) > on the support of ( for all v 
with \v\ < 1. 

We split the integral in (|15.10|) in two terms accordingly to (|15.11j) . Then the second term is 
real for e = 0. Now it remains to calculate the imaginary part of h(iu + 0), where 

f k 2 \p(k)\ 2 dk 
h(iuj + s) := / C(fc) ' , 7- (15.12) 

One can rewrite ()15.12j) as the iterated integral: over the surfaces T£ = {k e IR 3 : a(k) - b(k) = 
a, |.Do(&)| < 5} an d over a. Then we get 

k 2 \p(k)\ 2 dS 



and therefore 



h(iuj + e) = I U Q da, u(a) = / C(^) 3 !„ ; — 

' a + ie ' K J y r « v J 2a\V{a-b)\ 



, k 2 \p(k)\ 2 dS 



2a\V(a-b)\ 



This implies ffTpHJ) . since D Q (k) = a 2 - b 2 and |V-D (fc)| = 2a|V(a - b)\ on T w . U 
This completes the proofs of the Lemma 15.3 and the Proposition 15.1. 

Corollary 15.4 Proposition \l 5. 1\ implies that the matrix M~ l {iuj) is smooth with respect to uj 6 IR 
outside the three points uj = 0, ±/i. 

Remark 15.5 The proof of the Lemma fl5.3l is the unique point in the paper where the Wiener 
condition is indispensable. In Lemma 115.21 we use only that the coupling function p(x) is not 
identically zero. 
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16 Singular Spectral Points 



Recall that the formula ()13.21j) expresses the Fourier-Laplace transforms Q(X), P(X). Hence, the 
vector components of the solution are given by the Fourier integral 

{m) = hI^ M ' l ^{t) du (161) 

which converges in the sense of distributions. It remains to prove the continuity and decay of the 
vector components. The Corollary 115 .41 bv itself is insufficient to prove the convergence and decay 
of the integral. Namely, we need an additional information about the regularity of the matrix 
M- l {iu) at the singular points u = 0, ±/i and about some bounds at \u\ — > oo. We shall study 
the points separately. 

I. Consider first the points ±/i. 

Lemma 16.1 The matrix M _1 (zu;) admits the following Puiseux expansion in a neighborhood of 
rkji: there exists an e± > s.t. 

oo 

M' 1 ^) = J^Rki^T^ 2 , |w=F/i|<e±, weE. (16.2) 

Proof It suffices to prove a similar expansion for M(icu). Then ()16.2j) holds for M~ l {iu) as well, 
since the matrices M(±ip) are invertible. The asymptotics for M(iu) holds by the convolution 
representation in ()13.19j) : 

Hij(X) = (g x * djp, dip) (16.3) 
since g\ admits the corresponding Puiseux expansions by formula ([13. 14)1 . ■ 

II. Second, we study the asymptotic behavior of M~ 1 (X) at infinity. Let us recall that M _1 (A) 
was originally defined for Re A > 0, but it admits a meromorphic continuation to the Riemann 
surface of the function a/A 2 + p 2 (see Lemma fl3.3j) . 

The following proposition is a very particular case of a general fundamental theorem about the 
bound for the truncated resolvent on the continuous spectrum. The bound plays a crucial role in 
the study of the long-time asymptotics of general linear hyperbolic PDEs, |4Uj . 

Proposition 16.2 We can find a matrix R Q and a matrix-function R\(oS) such that 

R 

M- 1 Uuj) = — +R x {u) 1 b|>/z + l, io G 1R, 

us 

where 

\dlR x (u)\ < r 4, \u\>p + l, cugIR (16.4) 

for every k — 0, 1, 2, 

Proof By the structure ()15.1|) of the matrix M{iuj) it suffices to prove the following estimate for 
the elements of the matrix H{iu) := H(iu + 0): 

C 

IfifrH^iu;)! < T 4, LoeM, \u\>fi+l, j = 1,2,3. (16.5) 
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Let us rewrite (j!6.3|) as 



Hij(X) = {D-\\)d jP) d iP ), Re A > 0, (16.6) 

where -D(A) is the operator (jl3.4J) . and D~ 1 (X) is a bounded operator on L 2 (IR 3 ). Let us denote 
by Br the ball {x G 1R 3 : \x\ < R}. Estimate (|16.5|) immediately follows from a more general 
bound 

II^D-^iw + O)/!^^^ ^^11/11^,, c^gIR, H>/*+l (16.7) 

\U3\ 

which holds for every R > and all functions f(y) G L\ := {f(y) G L 2 (JR 3 ) : supp/ C B R }. 
Namely, by (jl.9j) the asymptotics ()16.5|) follows from the bound ()16.7|) applied to the function 
f(y) = dj P (y) G L\ with R > R p . The bound (|16.7|) follows from a general estimate [3HJ Thm 3] 
(see also the bound (A.2')] [21, Thm 8.1], [39, Thm 3]). 

III. Finally, consider the point u = which is the most singular. This is an isolated pole of a 
finite degree by Lemma fl3.3[ and hence the Laurent expansion holds, 

n 

M-\iuj) = ^L k u~ k - 1 + h(u), \uj\<e , (16.8) 

where L k are 6x6 complex matrices, Eq > 0, and h(u) is an analytic matrix-valued function for 
complex u with \uj\ < Eq. 



17 Time Decay of the Vector Components 

Here we prove the decay ()6.19|) for the components Q(t) and P(t). 

Lemma 17.1 Let X G Z v $- Then Q(t), P(t) are continuous and the following bound holds, 

C(p,v, d ) 



\Q(t)\ + \P(t)\< 



\t\) 3/2 ' 



t > o. 



17.1 



Proof Expansions (|16.2j) . (|16.4j) . and (|16.8|) imply the convergence of the Fourier integral (|16.1j) 
in the sense of distributions to a continuous function of t > 0. Let us prove the decay (|17.1j) . We 
know that the linearized dynamics admits the secular solutions without decay, see ()6.24j) . The 
formulas ()3.3|) give the corresponding components Qs(t) and Ps(t) of the secular solutions, 



Qs(t) 
Ps(t) 











t + 





dviPv 



'17.21 



We claim that the symplectic orthogonality condition leads to (|17.1jl . Let us split the Fourier 
integral (|16.1j) into three terms by using the partition of unity (i(uj) + (2^) + (3(00) — 1, uj G H: 



Q(t) 
P(t) 



— j e^GM + C 2 (w) + Cs(^))M- 1 ^ + 0) 
h(t) + I 2 (t)+I 3 (t), 



Qo 
p' 



duo 



(17.3) 
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where the functions Cfc(^) £ C°°(IR) are supported by 



suppCi 


c 


{u 


g IR 


£o/2 < \u\ < /2 + 2} 


supp C 2 


C 


{u 


e IR 


M > A*+ 1} 


supp C 3 


C 


{u 


G IR 


M < ^o} 



;i7.4) 



Then 

i) The function G C°°(JR) decays like (1 + |t|)~ 3 / 2 by the Puiseux expansion ()16.2|) . 

ii) The function I 2 (t) G C[0, 00) decays faster than any power of t due to Proposition 116.21 

hi) Finally, the function Is(t) generally does not decay if n > in the Laurent expansion (J16.8)) . 
Namely, the contribution of the analytic function h(ui) is a smooth function of t G H, and decays 
faster than any power of t. On the other hand, the contribution of the Laurent series, 

( Q P$) ) := i / elU3tUu) £ Lk{u ~ tE M ' (17 - 5) 

is a polynomial function of t G IR of a degree < n, modulo smooth functions of t G IR decaying 
faster than any power of t. This follows by the Cauchy theorem applied to the integral (jl7.5j) if 
we change the integral over uj G [—£ /2, e /2], where (3(00) = 1, by the integral over the semicircle 
e l9 e /2, 9 G [vr,0]. Let us note that the formula (|17.2|) gives an example of polynomial function 
arising from (j!7.5|) . 

We must show that the symplectic orthogonality condition eliminates the polynomial functions. 
Our main difficulty is that we know nothing about the order n of the pole and about the Laurent 
coefficients Lk of the matrix M _1 (ic<j) at u = 0. Our crucial observation has the following form: 

a) The components (jl7.2j) of the secular solutions form a linear space Cs of dimension dim £5 = 6. 

b) The polynomial functions in (jl7.5j) belong to a linear space Cl of dimension dim Cl < 6 since 
(Qo,^)eIR 6 . 

c) Cs C Cl since any function (|17.2|) admits a representation of the form (117. 5|) . The validity of 
this representation follows from the fact that the secular solutions (J6.24|) can be reproduced by 
our calculations with the Laplace transform. 

Therefore, we can conclude that 

C L = C s . (17.6) 

Let us show that the secular solutions are forbidden since Xq G Z v ^, and hence the polynomial 
terms in (jl7.5|) vanish, which implies the decay (|17.1|) . 

First, the constructed vector components Q(t) and P(t) are continuous functions of t > 0. 
Hence, the corresponding field components \l/(t) and IT(t) can be constructed by solving the first 
two equations of (|6.23j) . where A\ is given by (|4.9j) with w = v = v(ti) (see (118. 1|) below). 
Proposition 118. ll i) in the next section implies that X{t) G C(1R, £). 

Second, the condition X G Z Vj p implies that the entire trajectory X(t) lies in Z Vt p. This 
follows from the invariance of the space Z v ^ under the generator A VyV (cf. Remark 16. 6|) . In other 
words, X{t) = P v X(t). 

On the other hand, identity (|17.6|) implies that X(t) can be corrected by a secular solution 
Xs(t) s.t. the corresponding components Q&{t) and -Pa(^) of the difference A(t) := X{t) — Xsif) 
decay at the rate (1 + |t|)~ 3/2 . Note that P v A(t) = P v X(t) = X(t) since P v X s {t) = 0. 
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Further, the difference A(t) £ C(JR, £) is a solution to the linearized equation ()6.23|) . Hence, 
the corresponding norms of the field components of A{t) also decay like (1 + |t|) -3 ^ 2 that follows 
from Proposition I18.ll ii). Therefore, || A.(t) ||— ^ < C(l + |t|)~ 3//2 , hence the components Q(t) and 
P(t), of X(t) = P v A(t) also decay like (1 + |t|)~ 3/2 . ■ 



18 Time Decay of Fields 

In Sections 12-17 we denote by X(t) the solution to the linearized equation ()6.23|) with a fixed 
initial condition X . Here we consider an arbitrary solution X(t) = (ty(-,t),Tl(-,t),Q(t), P(t)) of 
the linearized equation. We shall prove a proposition which can be applied to the solution X(t) 
from previous sections as well as to the solution A(t) above. Let us study the field part of the 
solution, F(t) = (^(-, t), n(-, t)), solving the first two equations from the system ()6.23|) . These 
two equations have the form 

Ht) = ( a V , l.rr) Hi) + ( ° n ,A . rr, V (18-1) 



A - m 2 v ■ V J w V Qtt) ' 
We shall assume that the vector components decay, 

IQWI<f^p, *>0. (18-2) 
Proposition 16.71 is reduced now to the following assertion. 

Proposition 18.1 i) Let Q{t) £ C([0, oo); H 3 ) ; and F £ T. Then equation 1118.1}) admits a 
unique solution F(t) £ C([0, oo);jF) with the initial condition F(0) = F . 

ii) If Fq £ Tp and if the decay \18.ty) holds, then the corresponding fields also decay uniformly 
with respect to v: 

for \v\ < v with any v £ (0, 1). 

Proof Step i) The statement i) follows from the Duhamel representation 



F(t) = W(t)F + 



t > 0, (18.4) 



where W(t) is the dynamical group of the modified Klein-Gordon equation 

Ht)=( a W, L .\ V(*)- (18-5) 



A — m 2 v ■ V 

The group W(t) can be expressed through the group Wo(t) of the standard Klein-Gordon equation 

•(«)=( A _° m2 o )m (18.6) 

Namely, the problem (|18.6|) corresponds to (|18.5|) . when v = 0, and it is easy to see that 

[W(t)F(0)](x) = [W o (t)F(0)](x + vt), x £ 1R 3 , tel. (18.7) 
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Denote by W(x — y, t) and Wq{% — y, t) the (distribution) integral matrix kernels of the operators 
W(i) and Wo(t) respectively. Then ()18.7|) implies that 

W(x — y,t) — Wq(x — y + vt, t), x,yeH 3 , teTR. (18.8) 

The identity ()18.7j) implies also the energy conservation law for the group W(t). Namely, for 
(*(•,*), n(-, t)) = W(t)F(0) we have 

/Dn ( «, ( ) - . ■ V»(*,.)l' + |V»(«,.)I' + m'l*<«, = const, ( e R 

In particular, this gives that 

\\W(t)F \\F<C(v)\\F \\r, teM. (18.9) 
This estimate and ()18.4|) imply the statement i). 

Step ii) The statement ii) follows from the Duhamel representation ()18.4|) and the next lemma. 
Lemma 18.2 For any (3 > 3/2, v < 1 and Fq e Fp, the following decay holds, 

\\W(t)F \\^< H+^ Foh, t>0, (18.10) 

for the dynamical group W{t) corresponding to the modified Klein-Gordon equation MS. 0^1 with 
\v\ < v. 

Proof The lemma can be proved by general methods of Jensen and Kato [21] relying on the 
fundamental Agmon's estimate the bound (A. 2')]. We give an independent short proof for the 
convenience of the reader. 

Step i) The matrix kernel Wq(x — y,t) of the group W (t) can be written explicitly since the 
solution to (|18.6|) has the form (see |2*2*j ) 

*(., t) = [^R(t) * + R(t) * n ], n(-, t) = *(., t). (i8.il) 

Here R(t) = R(-, t) = R (-, t) + R m (-, t), and 



j , . 8(t — \x\) „ . . m JtimJt 2 — \x 



47rt ' mv ' ' 4vr y/t 2 -\x\ 2 ' 
where 

and J i is the Bessel function of order 1. From here and well known asymptotics of the Bessel 
function it follows that 

w cm) = o, \z\ > t, 

\d?W (z,t)\ < C(5)(l + t)- 3 / 2 , \z\<(l-5)t, 
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for t > 1, \a\ < 1 and any 5 > 0. From the last two relations and (j!8.8|) it follows that, for any 
v < 1 and e = -^jp, the following estimates hold for the matrix kernel W(,z, t) of the group : 

W{z,t) = 0, |z| > (l + v)t, (18.12) 
|o£W(z,t)| < C(v)(l +t)- 3/2 , \z\<et, \a\<l. (18.13) 

Step ii) Let us fix an arbitrary t > 1, and split the initial function Fo in two terms, Fo = Fq t + Fg t 
such that 

II^IU + ll^ll/3<CPolU *>1, (18.14) 

and 

F^(x) = 0, \x\ > |, (18.15) 

F'^{x) = 0, |x| < j, (18.16) 

where e > is defined in (fTTTTSjl . The estimate for W(t)i^ t follows by IfnOty . fTTHT^) and $E2Q : 

\\W{t)F'l t \U < \\W{t)F'l t \\ r <C\\F^\\ r 

< C^WFgJpO- + ty?< C 2 ( £ )||F ||^(l +t)-P, t>l. (18.17) 

Step Hi) It remains to estimate W(t)FQ t . We split the operator W(t), for t > 1, in two terms: 

w(t) = (i-CMt) + CW(t), 

where £ is the operator of multiplication by the function ((\x\/t) such that ( = ((s) e Qj°(IR), 
C(s) = 1 for |s| < e/4, ((s) = for \s\ > e/2. Since 

|^C(NA)l<c> W<1> *>i, 

and 1 — £(|a;|/i) = for |x| < et/4, we have, for t > 1, 

||(1 - CWt)*itll-/J < C 3 (e)(l + t)^||(l - CW*)fy |* < C 4 (e)(l + 0^1 1^)^*1 

From here, ijTCTl and (118.141) it follows that 

||(1 - C)W{t)FZ )t \\-f> < C 5 (e)(l + t)^!!^!^ < C B (e)(l + t)^||F ||^ t > 1. (18.18) 

Step iv) Thus, in order to complete the proof of Lemma 118.21 it remains to receive a similar 
estimate for (W(t)FQ t . Let Xet/2 be the characteristic function of the ball \x\ < et/2. We will use 
the same notation for the operator of multiplication by this characteristic function. From ()18.15|) 
it follows that 

(W(t)F^ t = (W(t) X et/ 2 F^. 
The matrix kernel W'(x, y, t) of the operator C,W{t)xet/2 is equal to 

W'(x,y,t) = ((\x\/t)W{x - y,t) X et/2(y). 

Since ((\x\/t) = for |x| > et/2 and Xet/2{y) = for \y\ > et/2, the estimate ()18.13J) implies that 

\d%W'{x,y,t)\ < C(v)(l +t)- 3/2 , \a\<l, t>l. (18.19) 
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The norm of the operator (W(t)x £ t/2 '■ —>■ T-p is equivalent to the norm of the operator 

(1 + \x\yP(W(t) X et/2(l + \y\)~ P :T^T. 
The norm of the later operator does not exceed the sum in a, \a\ < 1 of the norms of operators 

+ \x\)-P(W(t) X et/2{i + \y\r p ] : l 2 (wl 3 ) © l 2 (wl 3 ) -> L 2 (M 3 ) © L 2 (M 3 ). (18.20) 

From ()18.19|) it follows that operators (jl8.20}) are Hilbert-Schmidt operators since (3 > 3/2, and 
their Hilbert-Schmidt norms do not exceed C(l + i)~ 3//2 . Hence 

\\(W(t)F^ t \\„ p < C(v)(l +t)- 3 / 2 \\F^ t \\ p < C 7 (v)(l +t )-y 2 \\F , t \\ p , t > 1. (18.21) 

The last estimate above is due to (118.141) . Finally, the estimates (I18.21|) . (118. 18|) and (118. 17|) imply 
(HEHUD. ■ 

A Appendix: Computing Symplectic Form 

Here we justify the formulas ()3.5|) - (j3.7|) for the matrix fl. For j, I = 1, 2, 3 it follows from ()3.3|) 
and (jSH that 

Q(Tj, n) = (djip v , diir v ) - (djir v , d t tp v ), (A.l) 

^( r i+3, ti +3 ) = (d Vj ip v , d Vl ir v ) - (d v .7T v , d Vl ip v ), (A. 2) 

and 

Q(rj, r l+3 ) = -(djip v , d Vl ir v ) + (dj7r v , d Vl ip v ) + e,- ■ d Vl p v . (A. 3) 
Let us transfer to the Fourier representation. Set 

$(k) := (2tt)- 3/2 J e ikx i){x)dx. (A.4) 

It is easy to compute that 



Further, differentiating, we obtain 

2(kv)kj 

k 2 + m 2 — (kv ) 2 ri " Vj v k 2 + m 2 — (kv^ 2 



_ ? 2(kv)k i - , k 2 + m 2 + (kv) 2 ~ . 
a Vj i\) v = jrr^^v, o Vj ti v = ikj— - Tj-r^ipv, J = 1, 2, 3, (A.6) 



and 



6 ' V ' 



{\-v 2 fi 2 

Then for j, I — 1, 2, 3 we see from (jA.l|) by the Parseval identity that 



n(Tj,Ti) = -2i / kjkiikv)^] 2 dk = 0, (A.7) 
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since the integrand is odd in k. Similarly, by (|A.2j) . 

f kMkv)(k 2 + rn 2 + (kvf)\^ 
Q ^ Tl ^ = " 4 V (k 2 + m 2 -(kv) 2 ) 2 = °- 

Finally, by (|OJ), 

tt(Tj,T l+3 ) = J dk\^ v \ 2 kjk 

This completes the proof of ()3.5j) - ()3.7|) . 



+ 



2{kv) 



k 2 + m 2 + (kvy 

k 2 + m 2 — (kv) 2 ' A; 2 + m 2 — (few) 



B Appendix: Positivity of the Matrix F 

Here we justify the inequality used above in the proof of Lemma [15.21 

1 1 2 



m 2 + k 2 — (\v\ki + to) 2 m 2 + k 2 — (\v\k\ — uj) 2 m 2 + k 2 — (\v\ki) 2 
under the conditions ()15.7|) : 



> 



\v\ < 1, < \u\ < fi = mVl - v 2 . (B.l) 

Let us denote M 2 := m 2 + k 2 , r± := \v\k\ ± uj, and r := |t> \k\. Then the inequality reads, after 
cancellation by 2M, 

112 112 . , 

+ T^—T ~ 77 Z + 777— + 777-7- " 777": > «■ ( B - 2 ) 



M-r+ M-r_ M-r M + r+ M + r_ M + r 
The sum of the first three terms in (IB. 21) can be written as 



1 1 2 2w 2 

+ 777- - 77 = 7777— V777 ,,,r ( B -3) 



iV-w JV + w N (N + u)(N -uj)N 

where N := M — r . It is easy to check that iV ± uj > and N > under the conditions (jB.l|) . 
Hence, the sum ()B.3|) is positive (or positive infinite). Similarly, the sum of the last three terms 
in (jB.2|) also is positive. 
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